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ON THE DYNAMIC STRENGTH OF RIGID-PLASTIC 
BEAMS UNDER BLAST LOADS* 


Mario G. Salvadori, ! M. ASCE and Paul Weidlinger ,2 M. ASCE 
(Proc. Paper 1389) 


SYNOPSIS 


Upper bounds to the dynamic strength of simply-supported beams acted 
upon by a uniformly distributed blast pressure are obtained under the as- 
sumption of rigid-plastic behavior in bending and/or in shear (development 
of “plastic hinges” and/or “plastic slides”). 

The time variation of the blast pressure is approximated by a negative 
exponential time function with peak value at t = 0. Maximum allowable values 
of the peak pressure are obtained for beams of rectangular and I section. 
The results of the present analysis are valid for other types of impulsive 
forces. 


INTRODUCTION 


The dynamic strength of an elasto-plastic system under an impulsive load 
producing fairly large plastic deformations may be evaluated in a sufficiently 
accurate, but very simple — by assuming that the system is made of an 
ideal rigid-plastic material. 3) The solution thus obtained neglects the 
elastic part of the deformation and assumes that the work done by the loads 
is all expended in plastic flow. The response predicted by this type of so- 
lution is accurate only if plastic deformations are large compared to elastic 


Note: Discussion open until March 1, 1958. Paper 1389 is part of the copyrighted 
Journal! of the Engineering Mechanics Division of the American Society of Civil 
Engineers, Vol. 83, No. EM 4, October, 1957. 
a. Presented at a meeting of the American Society of Civil Engineers in 
Pittsburgh, Pa., October, 1956. 

1. Prof. of Civ. Eng., Columbia Univ., New York, N. Y. and Associate, Paul 
Weidlinger, Cons. Eng., New York, N. Y. 

Cons. Engr., New York, N.Y. 
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deformations. (4) But, the rigid, plastic approach may be used to determine 
upper bounds of allowable blast loads, even if the elastic deformations are 
not negligible, by introducing a criterion of strength based on allowable 
plastic strains, provided the time of maximum plastic displacement be large 
in comparison with the fundamental periods of the beam in bending and shear. 
The mathematical solution of blast problems is also greatly simplified by 
approximating the positive phase of the blast pressure between zero and in- 
finity by a single function, rather than by two or more functions requiring 
successive integrations of the differential equations. This can be done by as- 
suming that the uniform pressure p acting on the beam has a negative-ex- 
ponential time history: 


~~ 
p(T) = pe 


(1) 


where 7 is the time, p, the peak value of the pressure, and 7 the so-called 
decay constant. 

In order to compare pressures with different peak values and decay 
constants, the pressures of Eq. (1) will be assumed to have the same total 
impulse k: 


a 


at =k. (2) 


Whenever the pressure decay is rapid, the beam response is known to depend 
essentially on the impulse but not on the time variation of the pressure. 
Hence, the results of the present analysis, and in particular the graphs for 
the evaluation of maximum deflections and of maximum allowable pressures, 
may be used in connection with other pulse forms having the same total im- 
pulse. 

The beams considered are assumed to undergo plastic deformations: a) 
in bending, by developing “plastic hinges” at mid-span; b) in shear, by de- 
veloping “plastic slides” at the supports; c) in bending and shear, by develop- 
ing both plastic hinges and plastic slides. The beams are supposed to be 
braced laterally so that no lateral or local buckling may occur. 


The Development of Plastic Hinges 


A simply supported beam of length L, depth h. width b, and cross-section 
area A (Fig. 1) is loaded at time T = 0 with the uniform pressure of Eq. (1). 
We indicate by s, the uniaxial dynamic yield stress of its material (Fig. 2) 
and by = kj Ah its plastic section modulus. The plastic form factor k; 
varies between 1/4 (for a rectangular cross-section) and 1/2 (for an ideal I 
section). Under the assumption of rigid-plastic behavior in bending, the 
beam will not move unless the moment at mid-span,which takes its maximum 
value: 


1 2 
(2) 
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Simply supported beam under impulsive load. 
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Rigid-plastic stress-strain diagram. 
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at? = 0, is greater than the capacity plastic moment: 


Hence, motion will only start if pp is larger than: 


A 


If M, is slightly greater than Mp, the middle section of the beam is entire- 
ly plastic and the two halves of the beam rotate as rigid bodies about the 
supports, while the middle section moment has the constant value M,.. In other 
words, a plastic hinge develops at the middle section and the beam becomes a 
mechanism. 

For values of M, much greater than Mp an entire central portion of the 
beam may become plastic. This second possibility is considered in Sec. 5. 


Single-Hinge Displacements 


Referring the left half of the beam to coordinate axes as in Fig. 3, its angu- 
lar acceleration 6 is given by: 


Pp,» 
oo ~ (5) 


where m is the mass of the beam per unit of length and 


t= 


(6) 


is the non-dimensional time. Integrating Eq. (5) with rest conditions, 
= = O, and letting: 


r= (<4), (7) 


the center deflection of the beam, yc = Lg, is given non-dimensionally by: 
2 


+ t-1) - 


(8) 


yc becomes maximum, and the beam motion stops (y¥ c= O), at a time tp de- 
fined by the equation: 


1-eB-rt, 20, (9) 
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Fig. 4 gives tp as a function of r: 


tp = f(r). 
By Eqs. (8) and (9) the maximum value.) of Yc becomes: 
2 t 
r(i+s5-)] - (11) 


By Eq. (10), the right-hand member of this equation is a function of r only and 
may be written as: 


where ¢(r)is the function plotted in Fig. 5. 


Single-Hinge Moments and Shears 


Taking into account the inertia forces-m x 8, and the applied load, pb, per 
unit of length, the dynamic reaction Rp, shears Vp and bending moments Mp 
during motion are given non-dimensionally by: 


L 


gre 


(13) 


- 2°) + (2 e™*)(a - 32) (14) 


where: 


z = 2x/L (16) 
t 


Figs. 6 and 7 give vp and mp versus z for 1 e © varying between 1 and 10. 


These curves may be interpreted either as ‘shear and bending moment dia- 
grams at t = O for r varying between 1 and 1, or as diagrams of a given 


beam, i.e., for a given value of r, ast rows from zero on. The diagrams 
corresponding to tet = 1 are the static diagrams (straight line shear and par- 


abolic moment), * 


The Spreading of the Central Hinge 


An investigation of the first and second derivatives of mp with respect to z 
shows that mp is maximum at mid-span provided r > 1/3. Hence a single 
plastic hinge will be present at mid-span all through the motion for all 
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Non-dimensional shear in bending motion. 
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Non-dimensional moments in bending motion. 


q 
1.4 
| 
Le 
0.8 Y A 
/ 
| 
| | 
Y 
oO 
0.2 0.4 0.8 1.0 3 
| 
| 
j 


1389-10 October, 1957 
beams for which: 


8k_A 
= a 


and again equal to one at: 


(19) 


Inasmuch as the bending moment diagram was derived under the assumption 
that m, < 1 at all points (single hinge at z = 1), the present analysis breaks 
down theoretically when r < 1/3 and indicates that when this happens plasticity 
will spread to a central portion of the beam. 

However, an inspection of Fig. 8, (giving graphs of mp for 1/6<r<1/3), 
shows that, for 1/5<r<1/3, Mp remains practically horizontal between z = 1 
and z = zg, and that: 


* m,(2z,) 41.04 for r>1/5. (20) 


Hence, when 1/5<r<1/3, the single-hinge analysis gives a good approximation 
of the actual motion of a rigid-plastic beam of constant mass per unit length, 
whose central section is only 4% weaker in bending than the adjoining sections. 
In view of the uncertainties in the determination of the plastic flow stress s,,, 
this analysis may therefore be accepted as sufficiently accurate for 1/5<r<1. 


Final Displacements and Strains 


The assumption of rigid-plastic behavior allows the determination of dis- 
placements and rotations, and in particular of the final plastic angle at the 
central hinge (Fig. 9): 


f L 


(21) 


In order to evaluate strains it becomes necessary to take again into consider- 
ation the elastic behavior of the material. 

If the time of maximum deflection is large in comparison with the period of 
the fundamental bending mode of the beam, the influence of the strain rate on 
the magnitude of the final strains may be ignored and the final strains may be 
evaluated as if the dynamic bending moment diagram were applied statically 
to an elasto-plastic beam. This also implies that S, may be chosen equal to 
the static flow stress. 
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The stress-strain diagram for an elasto-plastic material is given in Fig. 
10, where € _ indicates the maximum elastic strain and S, the corresponding 
maximum elastic stress or flow stress. 


Indicating by M, the capacity elastic moment of the beam: 


Me = 5,8, (22) 
where S, is the elastic section modulus, the stress at the extreme fibers of 
all the sections for which m= M_ will be equal to S,- Hence the extreme fiber 


stress will be s, in a plastic zone from z = 1 (beam center) to z = z3, where 
Z3 is the only real root of the equation: 


Mp (z) = M,/M, (23) 
The ratio M-/M, = Se/Sp equals 2/3 for rectangular sections and one for 
ideal I sections. 
At the time tp of maximum displacement, the abscissa z3 limiting the 
plastic zone is defined by the equation: 


5 2(3-2) Bo 5 2(1-2) = M,/M, 


(24) 


Since zg decreases with r and increases with time, as seen by the curves of 
Fig. 7, the largest non-dimensional width of the plastic zone: 
dp =zl- 23 (25) 


is obtained for r = 1/5 at t = O, and may either be taken from the graphs of 
Fig. 8, or computed from the equation: 


2 2(3 - 2%) + (5) - = (26) 


dp is given in Fig. 11 as a function of M,/M, for r = 1/5. 
When r > 1/3 and motion occurs with a single hinge, the extreme fibers of 
the central section are displaced by an amount (Fig. 12): 


h h 
«Sand (27) 


Assuming that up, is evenly distributed over a length dp/L we obtain a lower 


bound for the final strain at the extreme fibers: 2 


eu 
fa, (28) 


When r< 1/3 the beam (Fig. 13) presents a kink at C and also a plastic 
curvature over the plasticized portion of the beam from C to D. But while 
the kink angle ©, varies substantially with the value of r and can only be de- 
termined by a complete solution of the rigid-plastic problem, the final angle 
6, between the rigid portions of the beam and the horizontal is a to 
the values of r, particularly if the motion is of the impulsive type. 3,4) This 
means that while the displacement u, is theoretically concentrated at the mid- 
span section in the single-hinge analysis, it is distributed over a length 4 if 
r< 1/3and a plastic zone is developed. Hence an estimate of & ¢ may still 
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be obtained by Eq. (28) for r<1/3, and this estimate will be a bound lower 
than forr> 1/3, since up is smaller for asingle-hinge than for a partially 
plastic beam, 


A Criterion for Allowable Peak Pressure in Bending 


A beam will be assumed to have reached its maximum allowable load when 
the final strain € , of Eq. (28) equals an allowable strain€ |: 


2u 
La, (29) 


The value of p, defined by Eq. (29) is an upper bound of p, since & ¢ is a lower 
bound for the strain. 

Substituting in Eq. (29), u, from Eq. (27), from Eq. (12), and (3 = p,/k 
from Eq. (2), we obtain a transcendental equation for the maximum allowable 
peak pressure pop in the form: 


B B 
Pop’ So Pop’, (30) 


a ™ 


and @ is the function graphed in Fig. 5. The constants B, and Bg are typical 
of the beam since the impulse k is assumed constant. Once a beam is chosen, 
i.e., By and Bg are given, Eq. (30) is easily solved by means of the graph of 
Fig. 14, giving XQ(X), where X = By/(op/80)- 


The Development of Plastic Slides 


The beam of the previous sections develops a maximum static shear at the 
ends given by 


L (32) 


Indicating by 7 , the dynamic yield stress in shear of the material, which by 
the von Mises criterion equals: 


8/3 (33) 


the plastic c2pacity shear of the beam is: 


‘5 A. (34) 
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Plastic longitudinal displacement at mid-span. 


FIG. 12 


Mid-span and total plastic angle. 
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Whenever the beam dimensions make it stronger in bending than in shear and 
Vo is larger than Vp, the end sections of the beam become entirely plastic in 
shear and the beam moves downward as a rigid body along two “plastic shear 
slides”, while the end shears maintain the value Vp: Plastic slides will thus 
occur if pp is larger than: 


“ote (35) 


Slide Displacements 


Referring the beam motion to the axes of Fig. 15, the beam acceleration is 
given by: 
(36) 


s 


s = V/V, (< 1). (38) 


A comparison of Eq. (37) with Eq. (8) shows that the law of motion of the 
center section of the beam in bending is identical with the law of motion in 
shear provided m be substituted for 2 m and r be equal to s. Hence the maxi- 


mum deflection in shear 5 8 is given by: 


and occurs at a time tg defined by the equation: 


1-eS-st #0. 


Q(s) appears in Fig. 5 and tg in Fig. 4. 


Slide Shears and Moments 


The shear and bending moment diagrams during slide motion are the static 
diagrams of a beam uniformly loaded, with end reactions equal to Vp- They 
are given in non-dimensional form by: 


P 
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Final Shear Strains 


When the time of maximum shear displacement is larger than the funda- 
mental shear period of the beam, shear strains may be approximated by as- 
suming that the beam behaves elasto-plastically in shear, as shown in Fig. 16, 
where v~ is the maximum elastic shear strain and € the corresponding flow 
stress. 

The portion of the beam ends in which the neutral plane fibers reach the 
strain?’ ., is determined non-dimensionally from the shear diagram of Eq. 
(41) by letting: 


(43) 


where V, is the maximum elastic shear. Since V,/V, = 2/3 for rectangular 
sections “and is more than 2/3 for I sections, the non- -dimensional width dg of 
the plastic shear zones varies between 1/3 and zero. 

For deep beams the elastic shear distribution across the depth of the beam 
differs from the distribution in shallow beams, and the maximum elastic 
strain does not occur at the neutral axis, as assumed above. But since the 
value of V, does not vary substantially with the depth of the beam, ifh<L/4, 
the evaluation of dg may still be obtained from Eq. (43). 


The value of the final meeeen’ fe then obtained by means of 5. and dg 
(Fig. 17): 
28 


at 


(44) 


A Criterion for Allowable Peak Pressure in Shear 


A beam will be assumed to have reached its maximum allowable load when 
its final shear strain 5 ¢ of Eq. (44) equals an allowable shear ween ii 


“et (45) 


s 


Substituting in Eq. (45) 5 ; from Eq. (39), and (3 from Eq. (2), we obtain the 


transcendental equation for the maximum allowable peak stress in shear in 
the form: 


Pog/T Pog! To (46) 
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5. = 2A 


= 


(47) 


The constants S; and Sp are typical of the beam. Equation (46) is easily 
solved by the graph of Fig. 14, where x = 8,/(P,,/7 ,) : 


The Possibility of Simultaneous Shear and Bending Displacements 


It was proved in Secs. 2 and 8 that bending motion can only start at a value 


of the pressure given by Eq. (4), and shear motion at a value of the pressure 
given by Eq. (35), but it is easy to see that, if the parameter: 


Pp 
oS ,min 1 L Ss 
h 


Pop,min 4/3 ky (48) 


which is typical if the beam, is less than 1, bending motion can never take 
place. 

In fact, the maximum dynamic moment during shear motion is independent 
of pp and has the value, given by Eq. (42) at z = 1: 


1 
Ms max. 


and this moment cannot be larger than Mp ifq < 1, since: 


(49) 


On the other hand, when q > 1 and Pp > PoB, min (and hence also py > 
Pos =PoB miaf®. shear and bending motions may occur simultaneously 
provided the maximum dynamic bending reaction Rp at % = O be larger than 
the shear capacity Vp, i.e., by Eq. (13) if: 


RL= 3M +M >VL, 
B P 


Solving this equation for M,: 


1 - sk. Ah 


and, substituting t | from Eq. (35) and s, from Eq. (4), the value of p, capable 
of inducing bending and shear motion is: 


min ~ Pop min * Pos ,min (4 - 3/a) . (50) 


Hence when q > 1, bending motion occurs for PoB, min< Po > PoS, min (4 
- 3/q), and simultaneous bending and shear motion develops for py > Pos, min 
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(4-3/q). Fig. 18 shows the various types of possible motion in a bending and 
shear beam. 


Simultaneous Slide and Hinge Displacements (q > 1) 


Indicating by Yp the acceleration of the left end of the beam and by % 


its angular acceleration, Fig. 19, the equations of motion in bending and shear 
become: 


from which: 


The acceleration of the center of the beam is Yo = Vp + B om or: 


2kM 
° 2 mL 

mL 


Integrating Eqs. (53) and (54) with rest conditions and letting: 


s, = (4 - 3/a)s r= (3 - 
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FIG. 19 


Beam in bending and shear motion. 
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and the relative displacement and speed of C with respect to A: 


“p> (¥, -¥,) = (q-1) t , (60) 
~ = (9-1) t . (61) 


The motion of A ends when ¥p = O, i.e., at a time ta defined by the equa- 
tion: 


-t 
= = O 
AA (62) 


and obtainable directly from Fig. 4 with s = sag. 
The absolute motion of C ends when Vo = 0, i.e., at a time tc defined by 
the equation: 


(63) 
and obtainable directly from Fig. 4 with r = rq. 
Since by Eq. (55) with q> 1: 
oh » (64) 
and t = f(r) is a decreasing function of r: 
ty < to (65) 
and A stops before C.* 5 
The maximum displacement © , of A is given by Eq. (39) with s = Sa: 
af (s,) 
A A’ * (66) 


The displacement S ¢ and velocity vc of C at t = ta by Eqs. (60) and (61) 


are: 
2 
2m 2 
3 (en) (67) 
2 af (68) 
Yo kr (q-1) t, 


*Since the translational motion stops before the rotational motion, there is no 
jump in momentum at t = t,. 
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Since at t = t, the motion of the beam ends stops, the additional motion of 
C for t >t, is single-hinge motion with initial conditions: 


‘ 
Under the assumption of small displacements, yc is added to § c to obtain 


the total displacement for C for t >t,. 
Letting: 


t' =t-t, (70) 


integration of Eq. (5), where 6 = 2y, /L with the initial conditions of Eqs. (69) 
gives: 


sen (e* -1+t') - 5 tt + le (4-1) (71) 
71 


2 mG -t 
ze A(l-e’) -rt' + 4r (q-1) t 


P 


(72) 


The motion of C stops when Ve = 0, i.e., at a time t{, defined by the equation: 


e*A (1 - e~*c) te = - kr(q-1) tas (73) 
Since ta is a function of r and q, t{, is also a function of these variables: 


= fo (r,a) ; (74) 


fc is plotted in Fig. 20 versus r with q as a parameter. 
Substitution of tt from Eq. (74) in Eq. (71) gives the maximum value} C 
of yc: 
t! 


2m -t 


= b 5 $,(r,a) 


4 
3 


(75) 


where @Q (r, q) is the function plotted in Fig. 21.versus r with q as a parameter, 
The maximum relative displacement d oof C with respect to A is the sum 
of 0 and 


P 


2 
= 


Y (r,q) = + 


a & 
| 
q 
2 
5 Ye 
| 
an 
or: 
= 
4 
| 
(76) 
where: 
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FIG. 20 


Stopping time of mid-span in bending and shear motion. 
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Maximum bending displacement in bending and shear motion. 
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is given in Fig. 22 versus r with q as a parameter. 


Slide-and-Hinge Shears and Moments 


The expressions for the accelerations of Eqs. (52), (53), and (54) allow the 
evaluation of the inertia forces during the slide-and-hinge motion and hence 
of the shears and the moments: 


(2) 2 BE ( 
= — = 4a (3z-1)(z-1) - 12 2(2-1) , 
BS (78) 


S 2 2 
m,(2) 27 = 2q z(1-z)~ + z~ (3-22) . (79) 


The shear and moment diagrams are given in Figs. 23 and 24,respectively. 

After the slide motion ends at t = ta, the shears and moments are given by 
the bending motion equations, i.e., Eqs. (14) and (15), in which tg +t' is sub- 
stituted for t. 


Spreading of the Central Hinge in Slide-and-Hinge Motion 


The single hinge analysis fails as soon as the moment diagram becomes 
minimum at z = 1; from Eq. (86) it is seen that this happens when q = 3/2. 
Since slide-and-hinge motion can only occur ifq > 1, this phase of motion only 
takes place for beams satisfying the limitations: 


. (80) 


When q> 3/2 the moment is maximum at: 


q-1 
and equals M, at: 


eq-e (82) 


The graph of Fig. 25 shows that if q ¢ 2 the moment is at most equal to 1.04 
Mp and that, hence, the single-hinge analysis may be considered practically 
valid, since it is rigorously correct for beams having a center section only 
4% weaker than the neighboring sections. 

After the slide motion stops, at t = ta, the single-hinge analysis theoreti- 
cally fails for: 


(83) 


As shown in Sec. 5, the analysis is still practically acceptable (within4¢g ac- 


curacy in the maximum moment) provided rz 1 etA, 
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Maximum displacement of mid-span relative to beam ends in bending and shear motion. 


FIG. 22 
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Non-dimensional maximum bending moment diagram in bending and shear motion. 
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It can be shown that additional slides cannot develop and the end slides 
cannot spread during bending-and-shear motion of beams with constant bend- 
ing strength under uniform pressure, although additional slides may appear 
(starting at the quarter points) in beams with weakened center sections which 
prevent the spreading of the central hinge. 


Criteria for Allowable Peak Pressures in Slide-and-Hinge Motion 


In the case of slide-and-hinge motion the peak allowable pressure may be 
determined either by shear or by bending. 

When shear governs, the non-dimensional width of the plastic zone, dg, 
is evaluated from the equation: 


Yas VL/M, (84) 


or from the graph of Fig. 23, and the final shear strain is equated to the al- 
lowable shear strain K a: 


= 
(85) 


Substituting dy from Eq. (66) and Ik from (2), the transcendental equation de- 
fining the peak pressure Pos in shear becomes: 


Pig/T 


oS ° 


' 2A . 
S} = (4-3/a) 
(87) 
Equation (86) is easily solved by the graph of Fig. 14, where X = S3/(P3,/> >) 


When bending governs, the non-dimensional width of the plastic zone Bis 
determined from the equation: 


Mps = (88) 


or from Fig. 24, and the final straine | is equated to the allowable strain E. 


(89) 
Substituting: 


2 
co) 
where: 
mt a 
wit 
a 
bel | 
u_ = q 
L c (90) 
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where d ¢ is given by Eq. (76) and /3 from Eq. (2), the transcendental equa- 
tion for the peak pressure pj, p in bending becomes: 


B B 
1 1 91 


where Bb is obtained from Eqs. (31) by substituting dj, for dg. Equation (91) 
is easily solved by the graphs of Fig. 26, where x+(X,q) is plotted versus X 
with q as a parameter. 

The lower of the two pressures pig and POR will be the allowable peak 
pressure in shear and bending motion. 


CONCLUSIONS 


The results of the present analysis allow the determination of upper bounds 
for allowable peak pressures of rigid-plastic beams under blast loads. Inas- 
much as the total blast energy is assumed absorbed in plastic flow and strains 
are supposed to be evenly distributed, within the limits for which this analysis 
is valid, actual elasto-plastic beams will not take higher pressures for the 
given allowable strains. 

The range of pressures, for which this simple analysis is valid, is limited 
because of the assumption of single-hinge development. 

The typical behavior of beams with respect to shear and bending motion is 
clearly brought out by this treatment. 

The analysis can be obviously extended to other types of beam supports and 
other types of impulsive loads. 
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Function for the evaluation of allowable peak pressure in bending and shear motion. 
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DEMONSTRATIONS OF PLASTIC BEHAVIOUR OF STEEL FRAMES 2 
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(Proc. Paper 1390) 


SYNOPSIS 


This paper presents the results of full scale tests which were made to il- 
lustrate the behaviour under load of structures designed by the plastic theory. 
The tests were made on beams, small scale frames, and a full scale pitched 
portal frame, all fabricated from standard rolled sections. Plastic analysis 
of moment redistribution was verified, and shear, instability and flange 
buckling effects were checked. 


INTRODUCTION 


The analytical basis for the Plastic Design of Steel Structures has been 
well established, and a large amount of testing has been carried out in con- 


nection with research.(1,2) Because of the unusual nature of Plastic Design, . 
there appears to be some need for full-scale demonstration tests on practical a 
structures and structural elements. However, the number of such tests fe 
which have been completed is still rather limited. The tests to be described é 


were performed to supplement a course of lectures on Plastic Design of Steel 
Structures given at the Royal Military College and Queen’s University in 
Kingston, Ontario, in May of 1956.(3) The tests were intended to demonstrate 
the characteristic behaviour under load of steel-framed structures and their 
elements - to show the accuracy of plastic methods of analysis and design, 
and to show the actual physical appearance of practical structures when 
loaded up to collapse. 


Note: Discussion open until March 1, 1958. Paper 1390 is part of the copyrighted 
Journal of the Engineering Mechanics Division of the American Society of Civil Engi- 
neers, Vol. 83, No. EM 4, October, 1957. 
a. Presented at a meeting of the American Society of Civil Engineers in 
Pittsburgh, Pa., October, 1956. 
*Now lecturer in Civ. Eng., The University of Glasgow, Scotland. 
**Asst. Prof. of Civ. Eng., Queen’s University, Canada. 
***Prof. and Head of Dept. of Civ. Eng., Royal Military College, Canada. 
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The tests are arranged to show: 


i) the moment-curvature relationship characteristic of structural 
steel members, 
ii) the fundamental idea of moment redistribution - the key to the plastic 
design of continuous structures, 
iii) the correspondence of predicted and observed collapse form and 
loads, 
iv) the effects of shear and instability, 
v) the behaviour of a practical rigid frame when loaded up to and past 
the working range, and on to failure. 


Testing Arrangements 


Most of the tests were performed on a test-bed made up of a pair of heavy 
rolled steel beams. Loads were applied by hydraulic jacks, and measured by 
pressure gauges located on the oil line near the pumps. The gauges were 
calibrated by jacking against the weighing system of a conventional testing 
machine, and the rate of testing was slow enough for the arrangement to be 
sufficiently accurate for the intended purposes. Some other tests were per- 
formed in conventional testing machines. Dynamometers with electric re- 
sistance strain gauges as sensing elements were used to weigh reactions and 
measure loads in a few cases. 

Slow rate coupon tests were made on each shape used and the results 
tabulated in Table 1 under the heading “Average Yield Stress.” These coupons 
were cut from 


a) a flange, extreme edge 
b) same flange, adjacent to the web, on the side of the web away from (a). 


The values of My, Mp and ultimate loads were calculated using the average of 


of the tensile yield stresses so found. 

Deflections and deformations were measured by means of fixed telescopes 
trained on scales fastened to the element undergoing test. No measurements 
of strain whatsoever were made on any of the test pieces. It may be noted 
that except for special matters such as fatigue, the actual stresses in a struc- 
ture are of no consequence providing the structure can support the required 
load without excessive deformation. It is an indication of the rational nature 
of plastic design that it is concerned with load-capacity and deformation of 
entire structures, rather than with stresses and strains. 


Moment Curvature Relationship 


The development of the plastic strength of steel frame structures depends 
on the formation of more or less fully plastic sections in the members. De- 
formation is concentrated at these sections, which behave as “stiff” hinges, 
and collapse occurs only when sufficient hinges have formed to transform the 
structure into a mechanism. In the theory the material is assumed to be 
ideally elasto-plastic, (i.e. perfectly elastic up to the yield stress and then 
able to strain indefinitely at constant stress), and it follows that the plastifi- 
cation of a standard structural shape will be as shown in Fig. 1. When the 
maximum fibre stress reaches the yield value it remains constant and, with 
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increasing strain, yield penetrates the section until full plasticity is 


reached (Fig. 1d). Curvature can then increase indefinitely at a constant 
moment. 


At any stage in the development of yielding the moment capacity can be 
found by taking moments of the stress blocks about the neutral axis. The 


curvature is easily obtained since, at the edge of the elastic core, the strain 
is just sufficient to cause yield, i.e., in Fig. Ic. 


de 
At full plasticity the fully plastic moment 
where Z Plastic Section modulus 


1.14 x Elastic modulus for WF and I sections 


Figure 2 shows the results of a test on a 1OWF 25 beam loaded at its third 
points. Curvature was measured by recording the change in distance between 
plates welded to the top and bottom flanges a distance of 10 inches apart as 
shown in Fig. 3. Section properties are given in Table 1. 

The early divergence of the experimental curve in Fig. 2 from the pre- 
dicted curve is due to local yielding caused by residual stresses in the beam. 
This was shown, during the test, by some flaking of mill scale at relatively 


low loads. In spite of this, the fully plastic moment capacity was realized, 
and the characteristic flat top of the M - @ curve demonstrates the ability of 
structural steel members to deform plastically. The extent of yielding is 
shown clearly in Fig. 3 by yield lines extending almost to the centre of the 
web. At the conclusion of the test the maximum strain in the extreme fibres 
was still below the work hardening level. Additional deformation would have 
finally caused strain hardening and a slight increase in moment capacity. 

It should be noted that this simple test also reveals some of the inconsis- 
tency of conventional design methods. At the load at which maximum fibre 
stresses are presumed to just reach 20 ksi, yielding has already occurred 
due to residual stresses. Thus, if the load capacity of such a member was 
in fact limited to the load at which maximum stresses actually reached 20 ksi, 
(as implied in the usual specifications), steel design would be seriously and 
unnecessarily restricted. On the other hand, the ultimate load as predicted 
by plastic theory is attained in spite of such factors as residual stresses and 
imperfections. 

For practical purposes it is most convenient to describe the moment- 
curvature relationship as shown in Fig. 4, using just the two tangent lines, 
and disregarding the curved transition between them. Although this approxi- 
mation is obviously satisfactory in defining moment capacity at full plasticity 
(large curveture), it is less accurate for determining deformations at inter- 
mediate stages of partial plasticity. 


Bending Moment Re-distribution 


For the realization of the fully plastic load capacity of an indeterminate 
structure, it is necessary that once the most highly stressed section reaches 
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full plasticity deformation must occur there (without loss of strength) thereby 
allowing other parts of the structure to take an increasing share of the load. 
The simplest demonstration of this behaviour is given by the propped canti- 
lever, which is only once indeterminate. 

A 12 ft. span, 8 B 13 beam was rigidly fixed at one end and simply 
supported at the other, and loaded at a point 40 inches from the fixed end. 
Lateral supports were provided at the load point and at a point 48 inches from 
the propped end. The reaction at the prop was measured by a resistance 
strain gauge dynamometer. The results of the test, curves of bending mo- 
ment and deflections, are given in Fig. 5. At a load of 28 kips, the support 
section reached full plasticity but continued to deform until, at a load of 31.5 
kips., full plasticity was approached at the load point. The structure was 
ther. transformed into a mechanism, and the deflection at the load point in- 
creased under constant load. 

The load at collapse can be readily calculated by plastic design methods. 
At collapse, fully yielded sections, i.e. plastic hinges, develop at A and B so 
that the bending moment diagram is defined, 


from Fig. 6b, Mp +? M, = We 


(L +b) 
ab 
For the section properties given in Table 1 and the dimensions given in Fig. 5 
W = 29,300 lb. It is seen from Fig. 5 that not until this load is reached do the 
deflections become excessive. For purposes of this calculation the plastic 
hinge at the fixed end A was assumed at a point one half inch away from the 
face of the support. This results in the segments a and b becoming 39-1/2 
inches and 103-1/2 inches respectively. Fig. 7 shows clearly the concentra- 
tion of deformation in the two “hinge” locations indicated in Fig. 6c. 

The combination of large shearing force and bending moment at the support 
section caused early yielding and a reduction in the fully plastic moment 
capacity at this section. Characteristic shear yield lines are shown in Fig. 8. 
On the other hand, and tending to compensate, the steep moment gradient at 
the clamped end leads to a very short yield zone there and large strains oc- 
cur. The resulting strain hardening caused the moment at A to continue to 
increase as the moment at B approached a constant value. 

In the test it was impossible to develop perfect fixity at the clamped end A. 
It may be noted that had A been perfectly fixed, Ma would have been equal to 
1.88 times Mpg in the elastic range. It may be seen from Fig. 5 that this dif- 
ference was not realized. Nevertheless, the fully plastic strength of the 
structure was developed. In general, the rotation and settlement of supports 
while perhaps affecting the order of formation of plastic hinges, have no ef- 
fect whatsoever in the ultimate plastic capacity so long as the required 
resistance at the support can be finally realized. 


and therefore W = M) 


Simple Collapse Mechanisms 


Rectangular portal frames, 8 ft in span and 4 ft high, fabricated throughout 
of 4B13 section, were tested to failure under three kinds of loading: 
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a) with a single vertical load, applied at midspan of the beam 
b) with a single horizontal load applied at the top of the columns 
c) with the vertical and horizontal loads combined. 


According to the Plastic Theory, failure of a structure requires that the 
structure as a whole or in part must be able to deform a finite amount under 
constant load. A mechanism with one degree of freedom is thus the minimum 
requirement for structural collapse. If a structure, or part, has “n” statical- 
ly indeterminate reactions, the occurrence of “n” plastic hinges forms a 
statically determinate system. The formation of one more plastic hinge pro- 
vides the necessary single degree of freedom. 

The plastic collapse modes associated with each of the three test loading 
conditions are shown in Fig. 9. The bending moment diagram associated with 
each collapse mode is characterised by a sufficient number of peaks of bend- 
ing moment of alternating sign, each equal to the fully plastic moment of the 
member, to provide enough plastic hinges to form a mechanism. The bending 
moment diagrams are shown in Fig. 10 and in each case the collapse load 
can be calculated from conditions of simple static equilibrium as shown in 
the figure. 

A pair of tests were carried out for each of the three loading arrangements, 
(6 tests in all), and the results are shown in Figs. 11, 12, and 13. A complete 
analysis of the combined loading case is set out in Appendix II; the same 
methods are also applicable to the simpler single load cases. 

In each test the behaviour of the frame was essentially elastic until the 
yield stress was exceeded at the most highly stressed section. Because of 
imperfect fixing of the column bases, the elastic stiffness of the frames (as 
revealed by the load-deflection curves) was somewhat less than elastic theory 
predictions. In fact, perfect rigidity of column footings is almost never 
found in practice. As a result, stresses predicted from elastic theory are 
often grossly inaccurate in such cases. It cannot be too often emphasized 
that such things as imperfect fixity, support settlements and residual stresses 
have no effect on the plastic strength of a structure. Following the load- 
deflection curves, it is seen that there is a marked change in deflection rate 
as each hinge is formed, and that between the formation of the penultimate 
and final hinges the curves become asymptotic to the ordinate of the predicted 
plastic collapse load. 

In the case of the frame carrying the single vertical load, the plastic 
hinges oecurred in the tops of the columns, where axial and bending effects 
are combined, rather than at the ends of the beam. There was no reduction 
in moment capacity due to this axial load, however. 

In the case of the test with the single horizontal load, hinges occurred at 
the bases of the columns and in the ends of the beam. 

Fig. 14 shows the frame subjected to combined loading after completion of 
the test. The flaked whitewash shows clearly the hinge locations and the 
carpenters square clamped to the “windward” knee shows the corner to have 

remained undeformed. This figure also shows the lateral bracing used in the 
tests. 

In each case, then, the methods of simple plastic design are seen to give 
an accurate accounting for behaviour and a realistic value for the collapse 
load. This must be contrasted with orthodox methods which do not indicate 
either ultimate strength or actual behaviour of a structure. 
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Effects of Shear, Local and Lateral Instability 


(a) Shear 


Although the simple plastic theory of structures is based upon the ideal 
elasto-plastic behaviour of a ductile member subjected to pure flexure 
(Figs. 2 and 4), members in actual structures are seldom found to be acting 
solely in flexure. The effect of shear on plastic behaviour is somewhat com- 
plex, but essentially the result is a reduction in plastic bending capacity due 
to yielding of the web under combined stresses. There is a marked diver- 
gence from linear elastic behaviour, and deflections increase rapidly as the 
flanges become fully plastic. 

The reduction in bending capacity due to shear is small except in cases 
where shears are relatively large. Methods for estimating the reduction 
have been developed by both the Cambridge and Lehigh research teams, (4,5) 
and the respective results are quoted below: 

(the symbols are defined in Appendix I). 


te Be fay | (1) 
2 7 


or in the form of an approximate design rule (2) becomes 


Mp 065 OW7 > (3) 


Fig. 15 shows the result of a test on an 8 ft-span 10WF25 beam loaded 
equally at points 21 in. distant from the end reactions. Marked on the load 
deflection curve are the capacities predicted by the above formulas. It is 
seen that the actual reduction in flexural strength due to shear is somewhat 
less than predicted, being only about 5 percent. Fig. 15 also shows that the 
combination of shear and flexure does not inhibit the plastic hinge rotation 
capacity. As full plasticity is approached, large deformations occur with no 
reduction in flexural capacity. Fig. 16, showing an end panel of the beam 
tested, reveals a combination of shearing and flexural deformation. 

It should be noted that as shear increases moment gradients increase 
proportionally. The steeper the moment gradient, the shorter must be the 
“length” of the plastic hinge and accordingly, the greater the curvature at the 
plastic hinge. It follows then that large shears are associated with sharp 
curvatures and a tendency to early strain-hardening. The result of this is 
that the reduction in capacity due to the shear is at least partly compensated 
for by the increased resistance due to strain-hardening. The theoretical 
results noted before do not take account of this phenomenon, and the result is 
as seen in Fig. 15 with the observed behaviour somewhat better than that 
predicted. 
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(b) Lateral Instability 


The problem of reduction of fully plastic moment capacity due to failure 
by lateral buckling has not yet been satisfactorily solved from the designers 
point of view. Work by both Horne(6) and the Lehigh Team(7) have produced 
values of the allowable unsupported length under a condition of pure bending. 
Plastic hinges usually occur at points where there is shear and a moment 
gradient, and the case of pure flexure seldom arises. The usual problem is 
thus that of a member mainly elastic but with localized fully plastic zones. 

A suitable solution for the lateral stability of such a member is not available. 

In practice, points of maximum moment, i.e. possible hinge locations, are 
braced laterally. Between these points, lateral bracing is provided according 
to the rules used in conventional elastic design. 

Fig. 17 shows the crippling lateral deformation of a 10 ft. span 10 WF 21 
beam which was loaded equally at two points 28 inches from the supports. No 
lateral bracing was provided, and the deformation seen in the photograph had 


resulted when the load was only 65 percent of the estimated plastic collapse 
load. 


(c) Local Instability 


The essential requirement of a “plastic hinge” is that large rotation must 
be realized at a nearly constant load. 

For a given total rotation, the required curvature (which is proportional to 
the maximum strain) depends on the length of the plastic hinge. With a steep 
moment gradient, a plastic hinge tends to be very short, and maximum strains 
may become quite large. Present knowledge indicates that if deformation up 
to strain-hardening can be tolerated, hinge rotations will be sufficient to per- 
mit full moment re-distribution and the realization of the fully plastic capaci- 
ty of practical structures. 

Beam sections commonly used - wide flange and I-shapes - sometimes 
tend to deform locally when strains become large. It is, of course, essential 
that the cross-section of a beam remains stable as a plastic hinge develops. 
To ensure a rotation capacity up to strain-hardening, it has been proposed(8) 
that sections be chosen such that 


flange breadth 17 


flange thickness 


Less rotational capacity is required of the web, and it has been proposed that 
for section stability up to the yield point 


web depth < 43 
web thickness : 


A 6 WF 15.5 beam of 10 ft. span was loaded at its third points and the full 
plastic moment of resistance attained without any signs of flange buckling. 
In this case the ratio of flange breadth to thickness was over 22. In cases of 
more localized deformation, for example, at the point of support of the 
propped cantilever (seen in Figs. 7 and 8) ultimate failure was accompanied 
by local buckling, though there was no reduction in strength. In some other 
tests as well, local strains were large and local deformations occurred, but 
not before collapse mechanisms had formed so that there was no reduction in 
carrying capacity. 
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Connections - Knee Joints 


Connections in conventional steel-skeleton construction generally serve 
only to transfer shear, a function which can be adequately fulfilled with 
riveted or bolted joints. In plastic design, however, more or less complete 
continuity of the structural frame is essential, and most connections must be 
able to transfer both shear and moment. Welding provides the best practical 
means of making such connections. In designing connections for “plastic” 
structures not only must strength requirements be met, but hinge rotation 
capacity must be assured in order that plastic redistribution of moments can 
proceed. Brittle or unstable behaviour must be avoided at all costs. If there 
is danger from either of these sources, the best solution is often to so 
strengthen the brittle or unstable element that it will function elastically 
while plastic deformations occur nearby. 

To date, the most common application of plastic design has been to single- 
storey sheds and frames. In these a rigid knee-joint is required. Although 
in some cases a built-up haunched knee may be suitable (in spite of its ten- 
dency to instability), the simple square knee-joint formed by welding directly 
together the leg and beam sections is most often used. Since the transfer of 
moment around the corner of the knee produces very large shearing stresses 
in the web, the web must be reinforced in most cases. Simple design rules 
have been proposed for this operation.(9,10) 

Fig. 18 shows the proportions of two knee joints which were tested, to- 
gether with the results of the tests. One knee joint, without web reinforce- 
ment, failed by shear in the web and the moment capacity realized was con- 
siderably less than that of the section used. The other was reinforced, and it 
is seen that its capacity was sensibly the same as that of the beam-section. 
It should also be noted that in both cases the load-deflection curves reveal 
satisfactory ductile behaviour, indicating that such connections are altogether 
suitable for use in structures designed on plastic theory. 


Behaviour of a Practical Frame 


Previous tests on full sized frames have been summarized by Beedle.(11) 
The following tests of two 20-foot span frames under side (wind) and vertical 
loading provide additional evidence of the applicability of Plastic Design 
methods. 

The general details of the frames are shown in Fig. 19; one feature dif- 
ferentiating these frames from previous tests lies in the use of practical 
erection joints (using high tensile bolts) in the rafters near points of low 
bending moment. As can be seen in Fig. 20 the joints did not reduce the 
structural stiffness of the frames. Lateral support was provided at the knees 
and at the vertical load points. 

The predicted collapse mechanism is that shown in Fig. 20b, with plastic 
hinges at B, C, E and F. The collapse load can be calculated by the principle 
of virtual work. 

Fig. 20 shows the measured vertical deflection at the ridge and the hori- 
zontal deflections at the knee. Using the section properties in Table 1 the 
predicted collapse load is: P = 16,250 lb. This was slightly exceeded. The 
results also show a divergence from linear elastic behaviour slightly below 
the load calculated to produce first yield at the lee knee. This might be 
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expected due to the effect of residual stresses and elastic stress concentra- 
tions. 

By the methods of Appendix II(b) an estimate can be made of the final de- 
flection at the point of failure; this is seen to give a reasonable, though con- 
servative, value. At the load in question, the frame is becoming seriously 
deformed and an accurate prediction of deflection is really not necessary. 

For the sake of comparison, the load at which the maximum stress 
reaches 26,600 psi. - the orthodox design combined live load - is compared 
in Fig. 20 with the design load found from plastic design methods. In the 
later case, the yield stress has been taken at its nominal value of 33,000 psi. 
and, for the combined wind and live loading, a load factor of 1.41 has been 
applied. It is seen that plastic design provides a clear, and safe, increase in 
possible working load of 40 percent. 

Fig. 21 shows the state of the frame at collapse, the flaking of the white- 
wash clearly showing the development of plastic hinges as predicted. After 
large deformation at the collapse load, final crippling occurred by local 
flange buckling at the knees as in Fig. 22 and by lateral distortion of the 
frames between points of lateral support. 

The test results show clearly that plastic theory accurately predicts be- 
haviour and load capacity. It is seen that plastic design methods are not only 
rational, but can lead to marked economy. 


CONCLUSION 


The tests which have been described present further evidence of the 
validity value of the plastic theory of structures. It is seen that the behaviour 
of practical structures - with the usual imperfections and residual effects - 


is accurately predicted by plastic theory up to and including the ultimate 
strength. 
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APPENDIX I 
Notation 

- Curvature 

E - Young’s Modulus 

Z - Plastic Section Modulus 

de - Half depth of Elastic Core 

Mp - Full Plastic Moment 

M’, - Fuli Plastic Moment, reduced due to shear. 

My, - Plastic Moment of Web 

= (yield stress x clear web area x clear web depth) — 4. 

T - Av. Web Shear 

Oy - Tensile Yield Stress 

a - Member Length, from point of max. moment to point of zero moment. 

dw clear web depth. 

d - overall depth of beam. 


APPENDIX II 


(a) Order of Hinge Formation 


The order of formation of plastic hinges and the loads at which individual 
plastic hinges occur, are of no importance in design when only the actual 
collapse load need be considered. An idea of the order may however be ob- 
tained by assuming a shape factor of unity. Thus 
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Moment at Yield = Fully Plastic Moment, or 
With this simplification, a complete analysis becomes fairly easy. 

The analysis shown in Fig. 23 has been carried out for the simple portal 
frame subjected to combined loading. Fig. 23(a) shows the elastic moment 
distribution. For My = Mp, the first plastic hinge will occur at E (the point 
of maximum moment) when, 


Increase of the loads by a factor “N,” will increase the bending moment at 
D to Mp, but at the same time the moment at E will exceed Mp. This obvi- 
ously cannot occur and the support moment must change by a value M! such 
that Mr remains equal to M). The bending moment is then as in Fig. 23 and 


31 
at D: 35 Wy L (Ny) + M! = M, 


33 
80 


and since = 


at E: W, L (Ny) - M! = M, 


N = 33/32 


Increase of the loads (W}) by a larger factor No will produce at C but 
Mg and Mp will both exceed M) so that the reactions must change by M! and 
H! to maintain all moments less than or equal to M). The new bending 


moment diagram is shown in Fig. 23(c) and i» 
at c: 24 w - HIL = 
31 
at D: Wy + M! - HIL = Mp a 
33 


Then No = 66/55. 


Finally, increase of the loads (W) by a factor Ng will produce a plastic 
hinge at A and reaction changes Ml, H1 and V! are necessary to maintain 
Mg, Mp and Mc equal to Mp. Ina similar manner N = 99/80. That is to 
say, the first hinge occurs at E when 


the second hinge occurs at D when 
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(a) elastic_bmd 
(Ist hinge) 


(d) collapse 
over de. 2mp+h'l= 


over cd: 2mp+v'l= wl(54) 
80 


ate: mp+m'= wl(33) 
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We = 1.03 Wy 
the third hinge occurs at C when 
Ws = 1.20 Wy 


collapse occurs when the fourth hinge forms at A 


= 1.24 Wy -32 


(b) Deflection 


It has been shown (4) that a reasonable estimate of the deflection when 
collapse is imminent can be obtained by applying slope-deflection equations 
to the statically determine structure which exists just before the formation 
of the last hinge. This analysis is carried out in Fig. 24 for the simple frame 
subjected to combined load. The slope deflection equation is most con- 
veniently used in the form 


= 9, 


4 L 
' 
where 6 2 end rotation of A due to lateral loads in the length L 
when both ends of length L are free to rotate 


A = relative deflection of the ends of length L, perpendicular to 
the original position. 


and all terms are positive clockwise on the member. 
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ANALYSIS OF CONTINUOUS BEAMS BY FOURIER SERIES 


Seng-Lip Lee,! A.M. ASCE 
(Proc. Paper 1399) 


SYNOPSIS 


This paper deals with the analysis of continuous beams by means of the 
expansion of arbitrary load function and the intermediate redundant reactions 
in infinite trigonometric series.2 The values of the redundant reactions cor- 
responding to any load distribution are determined by the application of 
Castigliano’s theorem. The same procedure is used to derive expressions 
for the influence lines for the reactions. 


Notations 
An, Aj, An; An Fourier coefficients for loads and reactions. 
a9, am Span ratio. 
Bh Fourier coefficient for deflection. 
Ch Defined by (10). 
Cj, Co, Cz, C4 Constants of integration. 
c Distance of concentrated load from left end. 
D4, Dy2, Dim --- Series defined by (17), (20) and (22). 
E Modulus of elasticity. 
F Defined by (16). 
I Moment of inertia. 
1 


Span length. 


Note: Discussion open until March 1, 1958. Paper 1399 is part of the copyrighted 
Journal of the Engineering Mechanics Division of the American Society of 
Civil Engineers, Vol. 83, No. EM 4, October, 1957. 


1. Asst. Prof. of Civ. Eng., Northwestern Technological Inst., Northwestern 
Univ., Evanston, Ill. 


2. “Mathematical Methods in Engineering” by T. v. Karman and M. A. Biot, 
McGraw-Hill Book Co., Inc., New York, 1940, p. 323. 
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Redundant moment on left end. 


Redundancy. 


Concentrated load. 
Arbitrary load function. 
Given load. 
Intermediate, or redundant, reactions. 
r Infinitely small length increment. 

Strain energy. 
Deflection. 
Independent variable, distance from left end. 
Deflection. 


Particular solution for y. 


dy, do, dm Displacements of intermediate supports. 


64 Angular displacement of the tangent to the elastic 
curve at left end. 


INTRODUCTION 


The deflection of a beam of constant cross section under arbitrary trans- 
verse load p(x) is given by the differential equation 


= px) (1) 


dx* 


where E is the modulus of elasticity and I the moment of inertia of the beam. 
The general solution of (1) is 


y(x) = C,+C,x+ C+ + ¥p(x) (2) 


where y, is a particular solution of (1). For a simply supported beam of 
span 1, the boundary conditions 


O 


yo) 
yu 
yo) 
yl) 


provide four equations for the determination of the constants of integration 
Ci, Co, C3 and C4. 


oO 


— 
| 3 
| 
My, ag 
m 
ped 
Yp(x) 
P | 
| 
= 
(3) 
= 0 a 
= 0 
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Elastic Curve 
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Let us consider the continuous beam shown in Fig. 1. It may be treated as 
a simply supported beam of span 1 subjected to the action of a given load q(x) 
and the reactions of the intermediate supports R,, Rg and Rm, which 
can be approximated by Fourier series of the type 


q(x) = A sin (4) 


< 
R, 2. An sin l 
met (5) 
se. 
Rin 2, in 
where 
2 
A, = Q(x) Sin dx (6) 
(al+r) 
=2 lim R Sin dx ~2Risin na, 
(a,l+r) 
A, = 2 lim Be sin MTX dx Sin (7) 
(a, l-r) 


(2,,1+ 1) 
Ar =2 Lim | Rm sin dx = 
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= > By sin (8) 


n=} l 


Substituting (4), (5) and (8) in (1) leads to 


(10) 


An, Al AL and _ are as given by (6) and (7), in which q(x) is positive 
if acting downward while R,, Ro, —— and Ry are positive if acting upward. 
Substituting (3) and (8) in (2) leads to the equation of the elastic curve 


sin 
(x) = in (11) 
(11) is defined if the values of Ry, Ro, —— Ry in terms of q(x) are known. 


Redundant Reactions 


Differentiating (11) twice with respect to x leads to 


2 oo 


y"(x) 2,74 C, Sin (12) 


n=! 


The strain energy in the beam is given by the expression 


l 5 


Application of Castigliano’s theorem results in 
wu _ 
dR, 


| 
= C (9) 
where 
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where d;, dy, —— and d,, are respectively the relative vertical displace- 
ments of Rj, Rg, —— and R,, with respect to the straight line passing through 
the points of end supports A and B, being positive if upward. Attention should 
be given to the difference between the sign convention just stated and that for 


the differential equation. Substituting (13) in (14), giving regards to (6), (7) 
and (10), leads to 


sin ntra, 


Sin nta, Sin ntra, 


sin nta, Sin nta,, 


+2 
sin’ 
Sin nta, Sin NTA, 
A, Sin ntra, 


A, Sin na, 


3 
4 
| 4 
(15) 
a 
me 
1 
where 
| 
= 
(16) 
co 
D, = 2, a 
| 
ik 2. | 
4 
a 
—_ 
is 2 n* 
ne} 


EM4 
= A,, Sin 
| 
| 
| 


D,, 
>, + Sin* | 
n=l 


Thus, Ry, Ro, —— and R,, can be determined by the solution of (15) for any 
loading conditions provided the relative displacements of the supports are 

known. Substituting the values of these reactions in (11) yields the expres- 
sion for the deflection of the beam. It should be observed that the series in 
(17) converge very rapidly as shown in the following example. 


Example 1 


Consider the two-span continuous beam shown in Fig. 2. Assuming that 
the center support does not yield under the load, d, = 0 and the first of (15) 


becomes 
R, D. (a) 


Substitution of known quantities in (6) and (17) yields 


A. Sin Sin 4 


Exact solution yields R, = 11/16P = .6875P. 
Substituting the value of R, in (11) leads to 


= (Sin .69 Sin sin 


| 
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Substituting these values in (a) leads to F 
R, = = = 690P ae 
n=! n* 2 8! 4 
3 
2P1 
p 
+ 


differentiating which we obtain 


cin NT) coe NTX 


The last two expressions provide a convenient means for determining the de- 
flection and slope of the elastic curve at any point along the continuous beam. 
It should be mentioned here that while y(x) and y'(x) converge rapidly, the 
convergence decreases with each differentiation. Knowing the redundant re- 
action, the bending moment and shear can easily be determined by conven- 
tional methods. 

Let us next consider the beam shown in Fig. 3, which may be treated as the 
limiting case of the beam shown in Fig. 1 as ajl approaches zero and Rj ap- 
proaches infinity while the product R,a,1 = M, remains constant. For this 
case, the first of (7) becomes 

 2nw 
A, = (18) 


and (15) becomes 


The first of (15) drops out in this case since all the coefficients become zero. 
Differentiating (13) with respect to My, giving regards to (18), leads to 


aU _ 
aM 
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(22) 


and 6, is the angular displacement of the tangent to the elastic curve at A, 


being positive if in the assumed direction of My. Under any loading conditions, 
Mj, Rog, —— and R,, can be determined by solving (19) and (21) simultaneous- 
ly. 


Example 2 


Consider the beam loaded as shown in Fig. 4. For this case 6, = 0, dg = 0 
and (19) and (21) become 


Ww sin dx = 20 - cos ) 


l n 3 


= (463 


co 
= 5A, 
| 
dake 
ac 
Tp’ 
(b) 
M—D.+R.D = | 
222-2 | he 
Substitution of known quantities in (6), (17), (20) and (22) yields . ; 
2 
2 
1.645 
= 1.015 W a5 
= 0.955W 


a 
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Substituting these values in (b) leads to 


517 Wi +0.4963R, - 5075 wl 
3,025 wl | 
Solution of (c) yields M, = .0233 wi? and Rg = .402 wi. a 


Influence Lines for Reactions 


The procedure just discussed can be used to derive expressions for the 
influence lines for the reactions as a unit load passes from one end of the 
beam to the other. To derive an expression for the influence line for Ry in 
the beam shown in Fig. 5(a), let us first consider the force system shown in 


Fig. 5(b). Displacing the point of support at Ry a unit distance, keeping 


yg = 0 and applying the reciprocal theorem to the two force systems lead to 


Py 
R, = =y (23) 


Thus, the elastic curve shown in Fig. 5(b) is the influence line for Ri For 
this case, (15) becomes? 


Solution of (24) yields 


D, D,, - (D,,)° (25) 


: 
DD, - (Dar | 


R' 


3. Note that * -1 and d, =Vo =0 


ibe 
Wed 
| 
| 
| 
4 
q 
or 4 
a —_ 
1 2 “22 = 
j 
= 
3 
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Substituting (25) in (11) and (23) leads to ; 

22 


Example 3 
In the beam shown in Fig. 5(a), if ay = i and ap = :, (17) yields 


Substituting these values in (26) leads to 


—(51sin” 45 sin ) sin (e) 


The ordinates of the wena” diagram for Rj, as given by (e) are tabulated 


in Table I for arguments of 7 , using only the first two terms of the expan- 
sion. 

Referring again to the beam shown in Fig. 5(a), an expression for the 
influence line for Ra shown in Fig. 5(c), may be derived in similar manner 
with slight modification. Observing that the relative displacements of the 
points of supports at Ry and Ry» shown in Fig. 5(c), with respect to the 


straight line passing through AB are respectively (1-a,) and (1-ag), it is 
evident that 


Ry = YOO = - = (1- (27) 


where w(x) is the equation of the elastic curve shown in Fig. 5(d). For the 
force system shown in Fig. 5(d), employing the same procedure discussed 
earlier results in the expression 


G)(x) = a.)D,,- )D,2} sin nwa, SI 
Substituting (28) in (27) yields 
{(! a,)D, -(i-a,)D }sin 
| 2 12 TUNA 


| 
) 
D, = Sin le 64 | 
— 1 nt... 2nn_ 3 ) 45 | 
| | 
| 2 2nt )= il 
ax, sin 3 4 64 ik 
aif 
aq 
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(29) is actually simpler than it appears to be. The next example illustrates 
this point. 


Example 4 


Consider again the beam discussed in Example 3. Substituting (d) in (29) 
leads to 


(14 sin -13 sin 207) sin 
(iasin 3 7 ) sin (f) 


The ordinates of the influence diagram for Ra, as given by (f), are tabulated 


in Table I for arguments of x/1, using only the first two terms of the 
expansion. 


CONCLUSIONS 


An advantage of this solution lies in the fact that (11), (26) and (29) are 
continuous over the intermediate supports. Consequently the deflection and 
slope of the elastic curve as well as the influence lines for the redundant re- 
actions are each given by one expression which takes account of all the spans 


The rapid convergence of the series used in the solution reduces the numeri- 
cal work to a minimum. 


TABLE I. Ordinates of Influence Diagrams for Rj and Ra 


| Se ASCE LEE 
4 
R = | x 
| 
| 
yi 
} 
a 
4+ 
| 0 1 1 1 2 5 1 
1 6 3 2 3 6 
Rj 0 -789 1.000 0 -.211 0 
Ra 1.000 .378 0 -.077 0 0448 0 1 

Loy 

¥ 
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GRAPHICAL SOLUTION OF EQUATIONS OF VIBRATIONS 


Wen-Hsiung Li,*A,M. ASCE 
(Proc. Paper 1412) 


SYNOPSIS 


A graphical method is presented for the analysis of dynamical systems of 
oscillations described by two simultaneous first-order ordinary differential 
equations. This method involves step-by-step graphical integration in a man- 
ner similar to Jacobsen’s general method of solving second-order ordinary 
differential equations. 

Examples of application are taken from hydraulic vibrations with non-linear 
damping force and arbitrary forcing function. Systems involving two sets of 
two simultaneous equations (two surge tanks connected to a conduit; and a dif- 
ferential surge tank) are included as examples. 


INTRODUCTION 


Many dynamical systems of oscillations are described by two simultaneous 
first-order ordinary differential equations.(1) The equations to be considered 
in this paper are 


dy 
F(x, t) (1a) 


and x = F(x, y, t) (1b) 


When analytical solutions are unattainable or difficult to obtain, approximate 
solutions can be obtained by stepwise integration. In this paper, a graphical 
method is presented for stepwise integration of Eqs. (1). 


Note: Discussion open until March 1, 1958. Paper 1412 is part of the copyrighted 
Journal of the Engineering Mechanics Division of the American Society 
of Civil Engineers, Vol 83, No. EM 4, October, 1957. 


*Associate Prof. of Civ. Eng., The Johns Hopkins Univ., Baltimore, Md. 
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When the functions fj and fg in Eqs. (1) are independent of the variable t, 
the equations can be reduced to 
dx 
which can be integrated by the Method of Isoclines.(1) The method presented 
in this paper can be used for equations where the functions fy and fy may 
involve the variable t. This method is similar to the graphical method of 
solving second-order ordinary differential equations advanced by Jacobsen.(2) 


A Graphical Method of Solving Two First-Order Equations 


Before considering two simultaneous first-order ordinary differential equa- 
tions in a more general form, first take the following linear equations with 
constant coefficients: 


constant M (3a) 


3 
constant N 


where m and n are non-zero constants of opposite signs. Eliminate the 
variable t from these equations: 


— dx (4) 2 

mx- M ny - N | 

Integrating, ‘ ‘ 
(x-x.) + = constant (5a) 


(5b) 


(5c) 


u 


(Sd) 


Here u is directly proportional to y for constant values of m and n. When 
Eqs. (3) are dimensionally homogeneous, u has the same dimensions as x. 

The solution of Eqs. (3) in the u-x plane (with u and x in the same scale) is, 
according to Eqs. (5), a circle with its center at (u,, x¢), as shown in Fig. 1. 
The radius of the circle is fixed by a pair of known values of u and x, such 

as (Up, X_) at tp. Any point on this circle represents the simultaneous values 
of u and x at a certain value of t. The value of dt required for the simul- 
taneous changes du and dx from one point to another at a distance ds apart 


| 
and dx 4 1 y = | lia 
4 
where 
ue = 
| 
é mM | | 
| 
| 
vag 
q 


ASCE LI 1412-3 


(see Fig. 1) can be shown to be related to the angle d@: With the angle 6 con- 
sidered positive in the counterclockwise direction in the u-x plane, 


de= 48 _(N-nu) dt _ (N- nu) dt _ 


Eq. (6) shows that the angle @ increases mA) with t for constant values of 
m and n. 


Next, consider two simultaneous first-order ordinary differential equations 
in a more general form: 


(7a) 


and 


ae 


(x,y, t) 


These equations can be rewritten as 


d 
+ mx = f,(x,yu,t) + mx 


M «x,y, t) (8a) 


d 
ny Noxy.t) 


The values of the constants m and n are chosen to be non-zero, one being 
positive and the other negative. For dynamical systems with positive linear 
restoring forces, this condition is automatically satisfied. Eqs. (8) differ 
from Eqs. (3) in that, instead of being two constants, M and N are functions 
of x, y, and t. Since x and y are functions of t alone, the values of M and 
N depend upon t. Approximate solutions for x and y can be obtained by 
integrating successively through short intervals At, during each of which M 
and N are taken as constants for approximation. Thus, according to Eqs. (5) 
and (6), this approximate solution consists of the end-points of a series of 
circular arcs (see Fig. 2). For each arc, the center is located by Eqs. (5c) 
and (5d), in which the values of M and N may be computed from the values 
of x, y, and t at the middle of the interval At. These values of x and y can 
readily be estimated from the graph by extrapolation. The angle A@ (see 
Fig. 2) for each arc, according to Eq. (6), is 


46 = -n/(- 2). at (9) 


Graphical Method of Solving Second-Order Equations 


F(x,y,t)+ ny 


The method advanced by J acobsen(2) for solving second-order ordinary 
differential equations can be deduced from the method presented above. Con- 


sider the equation dx 
= t) (10) 


ng 
AS 
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| 
4 ai 
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| 
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+px = O(x,y,t)+ = pi Fixy,t) (ta 
where 
* * (11b) 


and p is a constant. Comparing Eqs. (11) with Eqs. (8), we have in the u-x 
plane, according to Eqs. (5) and (9), 


Xe = F(x, (12b) 


ue = (12c) 


A8 = p. At (12d) 


The solution of Eq. (10) in the u-x plane, a phase-plane in this case, is given 
by a series of circular arcs, the centers of which are located on the x-axis. 
This is Jacobsen’s method of solving second-order ordinary differential 
equations. Examples of application of his method can be found in published 
papers.(3,4) 


Application to Hydraulic Vibrations 


While the graphical method presented in this paper is applicable to other 
dynamical systems described by Eqs. (1), problems from hydraulic vibrations 
are chosen for demonstrating the application of the method. The characteris- 
tics of hydraulic vibrations have been studied by others.(5,6,7,8) This graphi- 
cal method is used to find quantitative answers with given values of the para- 
meters involved. 

Three problems are chosen as examples: (1) The first problem is to find 
the maximum surge in a simple surge-tank as shown in Fig. 3a. (2) the 
second problem is to find the surges in two surge-tanks connected to a conduit 
as shown in Fig. 3b. This problem involves two sets of two simultaneous 
equations. (3) The third problem is to find the surge in a differential surge- 
tank as shown in Fig. 3c. In this problem, the motion is described by differ - 
ent sets of simultaneous equations under changing conditions. 

The equations used in these examples are those conventionally used by 
hydraulic engineers.(9) In these equations, water is assumed to be incompres- 
sible and the effects of water-hammer are thus excluded from consideration. 


Example 1. A Simple Surge-Tank 


The problem is to find the maximum surge of the water-surface in the 
simple surge-tank shown in Fig. 3a when the valve, originally closed, is sud- 
denly opened. The physical constants in ft.-sec. units are: 


| 
al 
| 
| 
4 
|. 
| 
and 
e 
| 
Al 
| 
Ae 
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Length of conduit, L = 27,000 ft. 
Cross-sectional area of conduit, A = 80 sq. ft. 
Cross-sectional area of tank, A; = 800 sq. ft. 
Reservoir surface above outlet, H = 550 ft. 
Head loss in conduit = c|v|v = 0.5| v| v 


When the valve is opened, the discharge at the outlet is given by 


Q = 40J/h, (for t>0) 


The equations governing the motion are 


h c|lvjv + 
and 
dh 15) 
These two equations can be reduced to the following form: 
dv, — 16 
clvlv) (16a) 
and 
dh_ AY 
dt At At (16b) 


Comparing Eqs. (16) with Eqs. (8), we have, with the variables h and v cor- 
responding to x and y respectively, m = g/L, n = -A/A¢, M = (H-.| v| v)g/L, 


and N = - Q/A;. Thus from Eqs. (5) and (9), we have, for each At in the u-h 
plane (u and h in the same scale), the solution given by a circular arc sub- 
tended by angle A@, where 


AL 
= vy = (17a) 


of center = H- cmMy = 


- @/AL (17¢) 


that is, 


Vv of center = = 0.5 Wh , (for t>0) (17d) 


and 


in radians At 


From past experience in stepwise integration for the maximum surge, except 
in cases where transient disturbances of short durations are involved, eight 


(17e) 


—— 
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4 a 
] a 
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4 
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or more steps usually prove satisfactory.(10) In this example, a value of 15 
seconds for At is used. Thus, A6@ = 9.38 degrees. 

In the u-h plane as shown in Fig. 4, construct the line (h = 500-0.5 |v | v) 
which according to Eq. (17b), locates the ordinate of the center of the arc. 
Also construct the line (v = 0.5Vh) which, according to Eq. (17d), locates 
the abscissa of the center of the arc. Thus, with an estimated mean values of 
v and h during an interval At, the center of the arc can easily be located. 
For example, to locate point 5 in Fig. 4 after point 4 has been obtained, 
extrapolate from the curve 0-4 to estimate the probable position of the mid- 
point 4' of the arc 4-5. Point I, the intersection of the (h = H - c|v |v) line 
and a vertical through point 4', gives the ordinate of the center of the arc 4-5. 
Point II, the intersection of a horizontal line through point 4' and the line 
(v =0.5 Vh), gives the abscissa of the center of the arc 4-5. Thus point III 
is the center of the arc 4-5. With a radius equal to 4-III and an angle A@ of 
9.38 degrees, point 5 is located. 

In constructing Fig. 4, no revision of the estimated position of the center 
was found necessary except for the first arc 0-1. 


Example 2. Two Surge-Tanks Connected to a Conduit 


Two simple surge-tanks connected to a conduit are shown in Fig. 3b. The 
problem is to find the surges in the tanks created by gradual closure of the 
valve at the outlet. The physical constants, in ft.-sec. units, are: 


Conduit lengths: Ly = 22,400 feet 
Lo = 22,200 feet 
Cross-sectional area of conduit, Ay = Ap = 80 sq. ft. 
Reservoir surface above outlet, H = 550 feet 
Initial steadv velocity in conduit = 10.0 ft. per sec. 
Cross-sectional areas of tanks: A; = 892 sq. ft. 
At2 = 715 sq. ft. 
Head loss in upper conduit = cq | vq | vy = 0.342|vq| vy 
Head loss in lower conduit = Cg | vg| v2 = 0.435| va| ve 
Initial surface elevation in upper tank = 550 - 0.342x102 
= 515.8 ft. 
Initial surface elevation in lower tank = 515.8 - 0.435 x 102 
= 472.3 ft. 


The closure of the valve at the outlet is completed in 90 seconds such that 
the discharge Q at the outlet is 


Q = 36.8(1-)/R, (for O<t<90) 


The equations governing the motion are 


(H-¢,lMly) 


(19b) 


@ 
: iad 
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Atz Atz 
By comparing Eqs. (19) with Eqs. (8), we have in the uy-h, plane (see ‘af 
Fig. 5), according to Eqs. (5) and (9): pe 
AL 
¥ = 7.9y, (21a) 


h, of center H- «lvly, (21b) 


u, of center = Aly 


A, 

that is, a 
V, of center = AzVe a 

A, (21d) 


AQ, in radians = At (21e) 


With At = 15 seconds, A6, = 9.76 degrees. 
Similarly, from Eqs. (20), we have in the u9-hg plane, 


Us = V2 = 8. 8 V2 (22a) 


h, of center 


Vz of center 


AQ, in radians = Ang At (22d) 


With At = 15 seconds, Ao = 10.9 degrees. 


3 
iy 
rey 
4 
aa 
; 
4 
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ws = H- 
Q 
= (22c) 
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In Fig. 5, the lines ( hy = H - cy| vy {v,) and ( hg = H - cg|vg|vg) are 
first constructed. The steps taken to determine points 10 are shown. The 
mid-point 9' of the arc 9-10 in each plane is first estimated by extrapolation. 
The center of the arc 9-10 in the uj-h, plane is ,ocated with point I on the 
(hy = H - cy|vy| vy) curve below (or above) point 9' and the value of v2 
taken from the ug-hg plane. The center of the arc 9-10 in the ug-hg plane 
is located by point III at a distance (H - hy) from point II which is on the 
(hg = H - cg|v2g| vg) curve below (or above) point 9'. The abscissa of this 
center is at vg = Q/Ag, which is zero for t greater than 90 seconds. 

The same problem was solved by the same graphical method with At = 30 
seconds instead of 15 seconds. Practically the same results were obtained. 
However, with the larger value of At, it became more difficult to estimate 
the position of the mid-points by extrapolation. Sometimes a revision of the 


estimated position was found necessary. No such revision was found neces- 
sary with At = 15 seconds. 


Example 3. A Differential Surge-Tank 


A differential surge-tank is shown in Fig. 3c. The problem is to find the 
surge in the tank due to a complete rejection of load on the turbines. The 
physical constants, in ft.-sec. units, are 


Headwater surface above tailwater, H = 442 ft. P 
Length of conduit, L = 43,200 ft. r 
Cross-sectional area of conduit, A = 251 sq. ft. ee 


Cross-sectional area of riser, A, = 201 sq. ft. 
Net cross-sectional area of tank, A; = 3,220 sq. ft. 
Top of riser above tailwater, H, = 500 ft. 

Port area, a = 40 sq. ft. 

Coefficient of discharge of port holes, C = 0.71 
Head loss in conduit = c| v| v = 0.382| v|v 

Initial water-surface elevation in tank and riser = 442 - 0.382x13.152 
= 376 ft. 


Upon complete load rejection at t = 0, the turbine gates will be closed in 10/3 
seconds such that the discharge Q through the turbines is given by 


.23 
Q= i7o[ 1- |r, (for O<t<2) (23) 


(In practice, Q may also be obtained from the turbine characteristic curves 
for given value of ht and gate opening at given value of t). 

The equations governing the motion may be divided into two groups: 
(1) When the riser is not filled to the top (hy < H,): 


Bh, = 2£(H- cy) (24a) 


— As. _ 
A, dt 
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LI 
A, = +Ca [t 2g(h,-h,) , (for (25) 


(2) When water is spilling over the top of riser (h, = H,): 


dy = 2(H -H,- clviv) (26a) 


(26b) 
Eqs. (26) should be used only when there is spilling over the top of the riser, 


> 


Immediately after the rejection of load on the turbines, the water-surface 
in the riser rises rapidly. Before the surface reaches the top of the riser, 
Eqs. (24) and (25) are applicable. In order to obtain a graph of proper pro- 
portion, it is desirable to rewrite Eq. (24a) as 


e+ = 2(H-ciiv+9h) 


By comparing Eqs. (28) and (24b) with Eqs. (8), we have in the u-h, plane, 
according to Eqs. (5) and (9), 


= 12:95v (29a) 


h- of center 0.| (H 24 clviv) + 0.9h, (29b) 
44.2 - 0.0382 |v|v + 0-9h, 


/ AL 
u of center + 


v of center 


| 
‘oo 
Bi; 
1a: x 
i.e. 
q 4 
| 
| 
- 
3 
(29c) 
= 0.908 J hr + A 
49, in radians = (29e) 
me | 
a4 Bins 
ke 


1412-10 EM 4 October, 1957 


In order to include the transient effect of closure, a value of 3 seconds for At 
is useg during the period when the water-surface in the riser is rising to the 
top. With At = 3 seconds, 46, = 16.6 degrees. The procedure of solution is 
shown in Figs. 6a and 6b. For example, to determine points 4 in the two 
planes, the mid-point 3' of the curve 3-4 in each plane is first estimated by 
extrapolation. The center of the arc 3-4 in the u-h, plane can be computed 
from Eqs. (29b) and (29d) with values of v, h,. and h; of the mid-points 3'. In 
the u-h; plane, Av is obtained from the u-h, plane, and Ah; is computed from 
Eq. (25): 


(30) 


ah, = + 0.0707 


This procedure 3 followed until h; reaches Hy, (at t = 21 seconds). 

When the water-surface in the riser reaches and remains essentially at 
the top of the riser, Eqs. (26) should be used. For the convenience of having 
a u-h¢ plane with the same u as defined by Eq. (29a) for the u-ht plane, Eq. 
(26a) should be rewritten as 


AG, 
+ 10 = (H -H, -clviv + 10-4 (30 


Comparing Eqs. (31) and (26b) with Eqs. (8), we have in the Ree plane, accord- 
ing to Eqs. (5) and (9), 


[AL 
lOA.g 


r 


= 12.95 v 


h, of center h, - (H, —H + clviv) (32b) 


h,- 0-6! Iviv - 92.8 


u of center = Ra (=O in this case)  (32c) 


AG, in radians = at (32d) 


With At = 20 seconds, AO = 6.9 degrees. The procedure of solution is shown 
in Fig. 6b. For example, to determine point 10, the mid-point 9' of the arc 
9-10 in the u-h; plane is estimated by extrapolation. The location of the 
center of the arc can be determined from Eq. (32b) with the values of v and 
hy of point 9'. This procedure is followed as long as Eq. (27) holds, i.e. 


Ah, > 0.0707 [500-h, -at 
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This required condition should be checked at each step. 

After point 11 (at t = 101 seconds), Eq. (33) is no longer satisfied. The 
water-surface in the riser is falling below the top of the riser. Eqs. (24), (25) 
and (29) should again be used. With At = 10 seconds, Eq. (29e) gives 

A6, = 55.3 degrees. The procedure of solution is the same as for the 
period immediately after load-rejection. 

The same problem has been solved by numerical integration by Rich (p. 152 
of Reference 10). His solution gives the same maximum elevation of water- 


surface in the tank (483.5 ft.) at the same time (262 seconds after load- 
rejection). 
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(a) Example 1. A Simple Surge-Tank 


Tay 


(b) Example 2. Two Surge-Tanks Connected to a Conduit 


(c) Example 3. A Differential Surge-Tank 


Figure 3. Hydraulic Systems in the Examples 
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Final position 
of equilibrium 


Min.h=444.7 ft. 
at t = 162 sec. ~— 


u in feet 
50 


v in ft./sec. (wv = u/9.17) 


Figure 4. Solution of Example 1. 
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for t=10 sec. 


=16.6° for At= 3 sec- 


100 


v im ft./sec. (v = u/12.%) 


Figure 6a. Solution of Example 3. (For u-h, plane, see Fig.6b) 
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10 


= w/12.95) 


1(t=101 sec.) 


> 
~ 
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u in feet 


Figure 6b. Solution of Example 3. (For unh, plane, see Fig.6a) 
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Journal of the 
ENGINEERING MECHANICS DIVISION 


Proceedings of the American Society of Civil Engineers 


THE STRENGTH OF VERY SLENDER BEAMS! 


E. F. Masur2 
(Proc. Paper 1413) 


SUMMARY 


The response of a slender beam to lateral loads and twisting couples is 
affected by the presence of bending moments in the plane of major stiffness 
much as the bending of beams may be influenced by the presence of axial 
forces. If, in addition, the major bending moments are statically indetermi- 
nate, and if the beam is sufficiently slender to admit relatively large lateral 
deformations, these may in turn affect the distribution of the principal bend- 
ing moments. The resulting nonlinear theory is the subject of this paper. 

After the establishment of the basic equations, it is shown that the inclu- 
sion of nonlinear terms in the strain-displacement relations corresponds 
generally to a stiffening of the structure as compared with the familiar linear 
theory. The (redistributed) major bending moments and reactions are shown 
to satisfy a minimum principle which represents an extension of the classical 
Castigliano Theorem. It is demonstrated further that, for increasing lateral 
loads and torsional moments, a limiting major bending moment distribution 
is approached asymptotically. For certain singular cases, the corresponding 
equilibrium configuration may not be unique, in which case the possibility of 
a snap-through (Durchschlag) phenomenon arises. 

The theory presented herein is corroborated experimentally with a fair 
degree of accuracy. Elastic behavior is assumed throughout. 


INTRODUCTION AND ESTABLISHMENT OF BASIC EQUATIONS 


In the present paper, a beam will be referred to as being “slender” when 
its moment of inertia ly about the (vertical) y-axis is much smaller than the 
moment of inertia I,; in addition, its torsion constant K will be assumed 


Note: Discussion open until March 1, 1958. Paper 1413 is part of the copyrighted 
Journal of the Engineering Mechanics Division of the American Society of 
Civil Engineers, Vol. 83, No. EM 4, October, 1957. 

1. This investigation was conducted under the sponsorship of the Office of 
Ordnance Research, U. S. Army. 

. Associate Prof. of Eng. Mechanics, Univ. of Michigan, Ann Arbor, Mich. 
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small in comparison with I,. An example of a slender beam is furnished by a 
rectangular beam, whose thickness t is small compared with its depth h. 

The response of such a beam to lateral loads and torsional couples 
in the presence of bending moments about the x-axis has been the sub- 
ject of numerous investigations.(1,2)3 It is shown in these papers how 
the bending moments influence, and often aggravate, the displacements 
of a slender beam. This is analogous to the behavior of beam- 
columns, whose response to loads in the presence of axial forces is 
well known. 

In all previous publications on the subject, the assumption is made more 
or less tacitly that the bending moments are either statically determinate or, 
in the event of statical indeterminacy in the major plane of stiffness, that they 
may be computed on the basis of the conventional linear theory. This may 
actually not be the case. In fact, if the lateral displacements u and the rota- 
tions § are sufficiently large, the introduction of nonlinear strain displace- 
ment relations may serve to modify the predicted moments. This question is 
explored in detail in what follows. It is shown that a redistribution of bending 
moments takes place, which serves to stiffen the structure relative to its 
predicted stiffness according to conventional theory. For example, if there 
are no vertical loads acting on the beam, the linear theory predicts, of course, 
vanishing bending moments everywhere in the absence of initial stresses. 
However, when certain nonlinear terms are included in the analysis, such 
moments do arise, and approach limiting values as the magnitude of the 
lateral loads and torques approaches infinity. 

Before proceeding to the analysis, it may be well to point out the limita- 
tions of the proposed theory. Actually, the nonlinearity is partial only in that 
the strain-displacement relations contain terms up to the second order, but 
ignore those of higher order. Physically this implies that the lateral dis- 
placements may be comparable to the thickness of the beam, but are still 
assumed to be small in relation to its length. Theories of an analogous na- 
ture are widely employed in connection with the analysis of structural ele- 
ments in which at least one dimension is much smaller than the remaining 
dimensions; the best known example is probably the plate theory of 
von Karman.(3) Since, furthermore, the effect of plasticity in this presenta- 
tion is ignored, it follows that the beam must be very slender to lend physical 
significance to the proposed theory. In the numerical example treated in a 
later section, appreciable deviations occur from linear theory at fiber 
stresses of about 30,000 psi in a beam whose depth-thickness ratio is 16:1. 

In what follows, let a beam of the type shown in Fig. 1 be subjected to 
lateral loads Aq and twisting couples At, in which q and t are given func- 
tions of z (measured along the axis of the beam), and A is a multiplier 
which is allowed to increase indefinitely. If Au(z) and A8(z) are the hori- 
zontal displacement and the rotation, respectively, u and 8 are governed by 
the equations of equilibrium. 


(Ely u")” + Pu" + = (1) 


(E-B")"- [(GK-P B']'+ (11 4P he) pap = t(Z) (2) 


3. Numbers in brackets refer to the Bibliography at the end of the paper. 


Reference (2) contains a fairly comprehensive list of publications on the 
subject. 
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where primes denote differentiation with respect to z. In these equations P 
represents the axial force (measured positive in compression), M(z) the bend- 
ing moment about the centroidal x-axis, and p(z) the given vertical load ap- 
plied at a distance a above the shear center S. E and G are the usual elastic 
constants, I” is the warping constant, p is the radius of gyration about S, and 
yo designates the position of S.4 Finally, k is a cross-sectional constant and 
is defined by 


k= - dA (3) 


which vanishes for sections symmetrical about the x-axis. 

Equations (1), (2), and (3) are the familiar equations of lateral bending and 
torsion,(4) although slight discrepancies are present.(5) Through a process 
which is entirely analogous to the one employed,(5) 5 it can be shown further 
that the vertical displacement v(z) satisfies the relationship 


V"=-(-17" + (4) 


in which the second term on the right side represents the effect of the non- 
linearity in the strain-displacement relations. M*(z) designates prestressing 
moments (if any) and has been incorporated for the sake of completeness. 

The inclusion of the nonlinear terms in Eq. (4) and the deletion of analogous 
terms from Eqs. (1) and (2) imply that the vertical displacements are much 
smaller than the horizontal and rotational displacements. This in turn re- 
quires that the beam be slender, as was outlined above. In the development 
of Eq. (4) it is assumed further that the force P be much smaller than the 
critical buckling force about the (strong) x-axis; in view of the slenderness 
of the beam, however, this represents no added restriction. 

The bending moments M(z) and M*(z) satisfy the equations of equilibrium 


M’= - pe) =o (5) 


and a set of appropriate natural boundary conditions. These, together with 
Eq. (5), determine M and M* uniquely (the latter trivially) if the structure is 


4. The assumption that S lies on the principal y axis introduces no signifi- 
cant loss of generality. 

5. The writer is indebted to Professor E. Reissner for suggesting, as an al- 
ternative, the derivation of these equations from large deflection plate 
theory. This process can be carried out for the technically most signifi- 
cant case of a thin rectangular beam and, with the exception of the intro- 
duction of the term (1-2)-1, leads to relationships which are identical 
with Eqs. (1) and (2) (properly simplified for this case, of course). How- 
ever, the form of Eq. (4) is slightly modified by this procedure. This 
discrepancy, which may affect the results somewhat in the presence of 

non-uniform bending, is due to a minor variation in the basic assumptions 

underlying plate theory as compared with the Euler-Bernoulli beam theory 
on which the present derivation is based. 
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statically determinate relative to its major bending moments. In that event, 
for given vertical load p(z) and axial force P, the “calibrated” functions u(z) 
and §(z) are also uniquely determined by Eq. (1) and (2) and the associated 
boundary conditions. Hence the total response Au and A increases in pro- 
portion to the total lateral load and torque Aq and At, respectively. 

A different picture is presented by a structure whose major bending mo- 
ments are statically redundant. In that case, the actual moment M(z) is dis- 
tinguished, among all moments satisfying Eq. (5) and boundary conditions, by 
being associated with a geometrically compatible vertical deflection v(z). 
Since the latter is related to M by means of Eq. (4), it becomes apparent that 
increasing lateral loads and torques, as expressed by increasing values of A, 
are accompanied by a redistribution of the principal bending moments. The 
equations governing this redistribution are developed in the remainder of this 
section. 


In a structure of nth degree of indeterminacy, the most general expression 
for M and M* is 


/1(Z) = + Ag Wag (Z) (6a) 
/9*(Z) = Ma (Z) (6b) 
in which m,(z) is governed, but not uniquely determined, by 


Ta 
Moe = -P(Z) ( 
and the same boundary conditions which apply to M. The set of self- 
equilibrated moments my(z) (r=1,2. . .n) satisfies 


Mp =O (7b) 


and equivalent homogeneous boundary conditions. The set of numbers A, 
(r =1,2. . .n), as yet unknown, will hereafter be referred to as “redundant 
parameters,” while Cr (r=1,2. . .n) constitutes a set of “prestressing para- 


meters.” Repeated Greek subscripts, as in Eq. (6), represent summation 
over the range 1 to n. It is finally convenient, and always possible, to select 
the set of functions m,(z) in such a way that the “orthonormality condition” 


/ (r=S8) 


EL ) (7c) 


is satisfied. 

The redundant parameters A; may now be determined by multiplying _ 
Eq. (7b) by v(z) and by integrating over the length of the structure; here v(z) 
represents any geometrically admissible vertical deflection function, that is, 
one which is sufficiently smooth and satisfies the geometric boundary and 
continuity conditions pertaining to the vertical deflection. In view of these 
restrictions, two integrations by parts lead to the relationships 


(8) 
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In particular, let v(z) be the actual deflection function v(z) satisfying 
Eq. (4). Then, by Eqs. (4), (6), and (7c), Eq. (8) is converted into 


[Me AZ 
Ix 

There is no generality lost in letting mo(z) be the actual bending moment, de - 

rived on the basis of linear theory, in the beam in the absence of prestressing 


and of lateral loads and torques. In other words, let the set of redundant para- 
meters A, vanish when all prestressing parameters Abd and the load para- 


meter A also vanish. When this is substituted in Eq. (9), it means that 


[ =O 


In view finally of Eq.(10), which constitutes an expression of the well-known 
and often used principle of virtual work, the “compatibility” conditions (9) 
become 


¥ 


Discussion of Large-Deflection Theory 


It is seen that there are as many Eqs. (11) as there are redundant para- 
meters. Since the bending moment M(z) is expressed in the form (6a), it is 
apparent that the nonlinear range is governed by the solution of the differen- 
tial equations (1) and (2), in which M(z) satisfies simultaneously the com- 
patibility equations derived above. The complexity of the resulting system of 
equations is therefore evident, and hence the necessity of solving it by inverse 
or trial-and-error methods. 

This is relatively easy for singly redundant beams. In that case, a value 
for At may be assumed arbitrarily at the start. With the bending moment 
M(z) so chosen, Eqs. (1) and (2) are solved and the solution (u,§) is inserted 
in (the single) Eq. (11), which in turn is solved for the load parameter A. 

If the value of A so obtained is real, and if the assumed moment M(z) is 
statically admissible in the sense defined later on in this section, then, by 
virtue of the uniqueness of the solution (proved also below), a point has been 
established in the load-response diagram. The latter may be completed by 
repeating this process for different values of Aj. 

For higher degrees of redundancy the process may become prohibitively 
laborious. However, some insight into the nonlinear behavior of the structure 
may be gained by means of a number of principles which are developed in 
what follows. In addition, it is possible that these principles may be instru- 
mental in the reduction of the numerical labor involved. 

To this end, let a quadratic form U be defined by 


2U(/4 B) = [cer +GKB'*) dz (12) 
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In Eq. (12) the first integral represents the bending and torsional strain 
energy, while the second integral can be shown to give the work done by the 
vertical load p(z) and the axial force P. Also, a bending moment M(z) will 
hereafter be referred to as being “statically admissible” if it satisfies the 
equations of equilibrium (5) and the stability condition 


(13) 


for all non-trivial functions (u, 8) which are geometrically admissible, that 
is, which are sufficiently smooth and satisfy the pertinent geometric boundary 
conditions. 


If further the “potential energy” V is defined by 


V(/4; 8) = U(M; B) 48)-wia4e) (4) 


in which the integral expression represents the work done by the lateral loads 
and the torsional moments, it can readily be shown that 


VIM, 4 A) VIM; B) (15) 


if M is statically admissible and if (u,8) satisfy Eqs. (1) and (2).6 It may be 
noteworthy that the equality sign in (15) implies trivial equality between (u, 8) 
and (u,8), provided only the inequality (13) is admitted. Conversely, if U=0 
for, say, (u1, 1), then the equality (15) applies to any (u,f) = (u,8) + c(uy, 4), 
in which c is an arbitrary number. The pair of functions (u;,A,) represents 
the fundamental buckling mode of the beam. " e 

With these definitions, let a bending moment M(z) = Mp+Aq@Mg be statical- 
ly admissible, relative to a given load, and let the functions (u,#) be associ- 
ated with M through the solution of Eqs. (1) and (2). Let further M and 
(u,8) represent, respectively, the correct bending moment and configuration 
for the same load; this implies, of course, that Eqs. (6a) and (11) are also 
satisfied. Since (u, 8) are certainly geometrically admissible, it follows that 


VOM; 4 B) (15a) 
similarly to (15). But 
VM; B)= 8) + (Aq -Ax) | -K 


= VIM; Au) (Aa - AK) 


in which the second equality is in consequence of Eq. (11). 

Moreover, if Eq. (1) is multiplied by u and Eq. (2) by £, and if several 
integrations by parts are performed, it is readily demonstrated on account of 
the boundary conditions that 


6. Since Eqs. (1) and (2) are the variational Euler equations of V, and 
since the actual natural boundary conditions are the variational boundary 
conditions of V, the correct solution (u,8) makes V stationary for all 
moments satisfying the equations of equilibrium. The additional minimum 
principle (15) represents an extension to the restricted class of moments 
which satisfy also the stability condition (13). 
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W(u, 8) = 2U(M; 6) (16) 


An identical relationship applies to the system (M; u,#). Thus, in view of the 
definition of V in Eq. (14), the inequality (15a) now becomes 


8) U4; th 8) (Xa Aa) An - (17) 
or: 


The last term on the right side of (17a) is positive definite. Hence, in the 
absence of initial bending moments (AF = 0), the following inequality is estab- 
lished: 


B)4 UM) S 455 MY) (18) 


with 
= 


Up represents the familiar expression for the strain energy in bending about 
the major axis. (18) follows from (17a) through the application of Eqs. (6a), 
(7c), and (10). An alternate statement of the inequality (18), which takes ac- 
count of Eq. (16), reads: 


8) + Us 4) (+4) (18a) 


These last two inequalities may be expressed in the following principle: 

Of all statically admissible bending moments, the actual one corresponds 
to the smallest value of A2U+U, (or, alternately, of A2W+2U,) if these values 
are determined on the basis of a deflected configuration which is related to 
the bending moments through the solution of Eqs. (1) and (2). 

This minimum principle, it may be well to emphasize, relates to the dis- 
tribution of the redundant moments rather than to the deflection configuration 
itself; that the latter is itself governed by minimum principles is well-known 
and was, in fact, utilized in the derivation of the present principle. It repre- 
sents an extension of the classical Castigliano “Theorem of Least Work” into 
the nonlinear range. Indeed, the Theorem of Castigliano represents a special 
case, which can be obtained from either (18) or (18a) by setting A equal to 
zero. 

To carry the discussion further, let the inequality (15a) be subtracted 
from (15). This results in the relationship 


7. Actually, the customary statement of the Theorem of Castigliano implies 
not a minimum, but only a stationary property of the strain energy in bend- 
ing. The present principle can be broadened in a similar fashion. In fact, 
if M satisfies the equation of equilibrium (5), but not necessarily the sta- 
bility condition (13), the quantities dealt with in (18) and (18a) can be shown 
to be stationary. In other words, if (u,8) as solved from Eqs. (1) and (2) 
are differentiable functions of A,y, then the satisfaction of Eq. (11) ensures 
that the derivative of (A2U+U,) or of (A2W42U,) with respect to A; 
vanishes. 
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If now both (M:u,§) and (M;u, 8) are assumed to satisfy the conditions of com- 
patibility (11), it follows that 


or, after some rearrangement of terms, 


Obviously, the inequality above is impossible for real values of the redundant 
parameters. On the other hand, the equality is possible only if M(z) and M(z) 
are trivially identical. This in turn implies the following principle: 

If there exists a statically admissible bending moment which satisfies the 
conditions of compatibility (11) relative to a configuration (u,8) which is a 
solution of Eqs. (1) and (2), that bending moment is the actual bending moment 
and is unique. 

In general, it follows from the uniqueness of the M(z) that u(z) and A(z) 
are also unique, as discussed before. However, singular cases may occur in 
which, owing to a condition of “neutral equilibrium,” two equilibrium con- 
figurations become possible. This was already hinted at following the estab- 
lishment of the inequality (15). A detailed discussion of this question is con- 
tained in the Appendix as well as in later sections in connection with the 
numerical example and the report on the experiments. 

Returning to the minimum principle established above, it is clear then that 
any statically admissible bending moment furnishes an upper bound to the two 
functions with which the principle deals. In what follows it will be the object 
to derive also a lower bound and to discuss a possible application of these 
bounds toward the simplification of the computational labor. a 

To this end, let a “kinematically admissible” bending moment M(z) be de- 
fined as one which satisfies the equation of equilibrium (5) (and which, there- 
fore, can be expressed in terms of Eq. (6a) subject to Eq. (7) ), which is com- 
patible in the sense of satisfying Eq. (11) for some geometrically consistent 
configuration (u,8), and which finally does not violate the inequality® 


U(M; 4, 2) <o (19) 


If this inequality is now subtracted from the inequality (15), in which the actual 
state (M;u,8) is compared with the state (M;u, B ), then, in view of Eq. (16), 


Ma ( B- WCU, @) =-weu, 6) (20) 


where jy; represents the set of redundant parameters designating M(z). The 


integrals in (20) may be replaced in line with Eq. (11); after some rparrange- 
ment of terms, this leads to 


8. The concept of “kinematic” admissibility follows from the fact that if th 
equality sign in (19) is satisfied the external work equals the “internal °» 

work.” Note, however, that the equations of equilibrium (1) and (2) need 

not be satisfied. 
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in which the last term, as before, is positive definite. Hence, by the same 
argument as the one employed previously 


or, in view of the relations (16) and (19) 


Since the actual bending moment M(z) is also kinematically admissible, it fol- 
lows then that: 

Of all kinematically admissible bending moments, the actual one corres- 
ponds to the smallest possible value of the functions a*w+2U, or 
NU + U,. 

It appears that this lower-bound principle may be useful in estimating the 
nonlinear response of the structure. In fact, let an upper and a lower bound 
be found for a given value of the load parameter A; an error estimate has then 
been established. Since the bounds can be narrowed down arbitrarily, a cor- 
rect solution may presumably be approached in this fashion. 


Asymptotic Behavior 


As the load parameter A grows beyond bounds, the bending moment M(z) 
and with it the deflection mode (u,f) usually approaches a limiting condition. 
This asymptotic behavior is discussed in this section in regard to both the 
governing equations and the appropriate energy principles. 

For given load p(z), it was shown that the stable domain (that is, the one in 
which the inequality (13) is satisfied) is generally a closed region in a space 
in which the coordinates of a point are given by the redundant parameters A, 
associated with M(z) through Eq. (6).(5) As p approaches an ultimate ig 
P,, this region shrinks to a point; for p>py, no statically admissible bending 
moments are available any more. In the present paper, p is assumed to be 
fixed and less than p,,. Because of the boundedness of the associated stable 
domain, it follows therefore that all the redundant parameters A, remain 
finite as the load parameter A goes to infinity. 

This determines the form of the governing equations of the limiting state. 
In fact, if A approaches infinity in Eq. (11) while the left side remains 


bounded, the asymptotic deflection mode (u \ B »? is seen to satisfy the set 
of equations 


together with the equations of equilibrium (1) and (2). It may be of some 
interest to note that this system of equations does not contain the prestressing 
parameters »r In other words, the “final” state is independent of whatever 


initial stresses may exist in the structure (due to settlement of support, tem- 


perature gradients, etc.), although the “history” of the structure does display 
such dependence. 
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By the same token, an asymptotic minimum principle can be derived. For 
if A is permitted to approach infinity in the inequality (18), it follows (for 
=0) that 


; %, Be) = %, Bes SU(M: 4) 


Wty, bo) Swed (23a) 


in which (M;u, 8) represents, as before, any statically adinissible bending mo- 
ment and associated deflection mode. The first equality in (23) follows from 
Eq. (22). Accordingly: 

The limiting bending moment is characterized by making U (or W) smaller 
than does any other statically admissible bending moment. 

In general, this minimum is in the interior of the region of stability. Itis 
approached asymptotically, and the corresponding configuration is unique. 

The existence of such an interior minimum is proved in the Appendix; also 
explored are certain singular cases, in which U may assume a minimum on 
the boundary of the stable domain. Suffice it to state here that in that event 
the limiting bending moment may be reached for a finite value of the load pa- 
rameter A. Also, the deflection mode need not be unique in that case, as 
would appear reasonable in view of the limitations on the uniqueness principle 
established in the preceding section. 

As before, the upper-bound principle expressed through (23) has a counter- 
part in the form of a limiting lower-bound principle. This is derived readily 
by considering (21), with A going to infinity. If then (u,, By) represents any 
pair of functions which, in addition to being otherwise acceptable, satisfies the 
set of equations (22), and if this “collapse mode” furthermore satisfies the 
condition9 


ZU 4, be) - Wh, &) =O (24) 


U (Me; Uc, Be) U (Mo; Ue, br) (25) 


W (te, Be) fo) (25a) 


In other words: 
Of all collapse modes, the actual one corresponds to the largest value of 
U (or W). 


Numerical Example 


In this section, the principles and equations of the preceding sections are 
applied to an illustrative example. For this purpose, consider a beam of 


9. Unless one of the terms in Eq. (24) vanishes, this condition can always be 
satisfied by selecting the amplitude and sign of the assumed collapse mode 
properly. This follows from the fact that U is quadratic, but W linear in 
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length L which is of thin rectangular cross section (I” = 0) and of second de- 
gree of redundancy relative to bending in its major plane by virtue of being 
elastically restrained at both ends; however, only one degree of indeterminacy 
need be considered here owing to the complete symmetry of the problem. 
Simple supports are provided at both ends, so far as lateral movement is 
concerned; a single lateral force of magnitude A is applied halfway between 
supports and at an eccentricity e above the center line. If, as before, the 
total response is denoted by (Au,Af), and if m is the restraint moment, then, 
in the absence of the vertical load p(z), Eqs. (1) and (2) and the boundary con- 
ditions take the following form: 


Ely mp" = 0 o< 
COSZE 4) 
Ulaj= ula) (26) 


(4) = (4 )- MB %) + G=o 


In Eq. (26) the second set of boundary conditions follows from the symmetry 
of the problem; hence the solution for (u,8) need only be considered for half 
the beam, that is, from the left support (z=0) to the center (z=L/2). 

For the case under consideration, the compatibility condition (11), in the 
absence of prestressing, becomes 


m=-2 (EZ) ¥ (27) 


where 


Cy 
represents the degree of end restraint and varies from zero (no restraint) to 


unity (full fixity). C denotes the “spring constant” and is measured in in.-lb 
per radian. 


For the range between 0 and L/2, the solution of Eq. (26) is given by 
A ga = (I- 4a) + 22 
Cos Z 


(28) 


in which 
A 4° é 
3 = 


where my, identifies the moment associated with neutral equilibrium, and ? 
is a dimensionless constant. 


If Eq. (28) is substituted in Eq. (27), this can be shown to lead to the rela- 
tionship 


= 
a 


|, 
tt j 
4 
q 
| 
C 
| 
. 


with 


This equation can be made dimensionless. In fact, let the dimensionless ratio 
Ww be defined by 


where AL my (30) 


with S, and S, designating, respectively, the section moduli about the y and 

x axes. It is readily apparent that Sy represents the maximum fiber stress 
in lateral bending, while 8; is the maximum fiber stress associated with the 
buckling moment m4. With these definitions, and in consideration of the re- 
lationship between the moments of inertia and section moduli for rectangular 
sections, Eq. (29) now takes the simple form: 


(8¢ w? = Yeca) (31) 


As w (and hence the lateral loadA) approach zero, both numerator and denomi- 
nator on the right side of Eq. (31) also approach zero. For very small values 

of w, Eq. (31) can therefore be replaced, through the usual limit procedure, 

by 

1?) w* = 3K Cx << (31a) 


Conversely, as the lateral load becomes increasingly large, it follows from 
Eq. (31) that 


/im 


F(x)=0 (31b) 


The best way to establish a functional relationship between @ and w is 
probably to assume a value for the former and to solve for the latter by 
means of Eq. (31), although, for small values of @, the inverse procedure may 
be followed through the use of Eq. (31a). The asymptotic magnitude of o (for 
very large lateral loads) may be obtained from Eq. (31b); it can readily be 
verified that the smallest positive root of Eq. (31b) is given by a =1/¢ for¢@ 
exceeding 77, while for ¢< 7, the smallest root must be computed through a 
trial-and-error process. 

With the w-a relationship thus established, the maximum lateral deflection, 


- which occurs at the midpoint of the beam, is found to be 
Umax = AC/ San % Ky (32) 


In Eq. (32), the quantity A was defined in Eq. (28), while xy represents the 
ratio between the computed maximum deflection and the equivalent value 
determined on the basis of the linear theory, which does not predict the 
development of restraining moments m (that is,@=0). For very small 
values of a, Eq. (32) can be approximated by 


7 
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Ky=/- =/- ) CA (32a) 


in which the second equality follows from Eq. (31a). 
Similary, the maximum rotation occurs also at the middle of the beam and 
is governed by the equation 


CB max - ten = ky (33) 
Ks = ( [% - (/- Oa) Pan % 


in which, similarly to the preceding equation, kg is defined as the ratio be- 
tween the computed maximum rotation and the corresponding rotation com- 
puted on the basis of the linear theory. As before, for very small values of 
a, the second equation of Eq. (33) can be approximated by 


Ke =/- (Yop) = (7g) (A<< (33a) 


Of some interest, finally, is the lateral deflection of the point of application 
of the force itself. This is found by adding Eqs. (32) and (33), or 


2 
max = A [ tan - (/-29a) =K (14/26?) 


(34) 
-/ 
K = (244) 3)(/412 $2) $4)? - % | 
Here, X is defined analogously to K,, and Kg. For small values ofa, this 
goes Over into 
2 
K =/- /-(%) 4128) (34a) 


It may be worth noting that, from Eq. (34), dx/da can be shown to be pro- 
portional to F(a). In other words, for very large values of w, dx/da vanishes. 
Since furthermore the second derivative is positive, it is seen that the limit- 
ing value of x is an absolute minimum if only stable values of @ (lying in the 
closed region between -7 and +7) are admitted in comparison. This, however, 
is not surprising. Since x is a measure of the total work done by the external 
force (all other terms being independent of a), this property of x is a natural 
concomitant of the minimum principle expressed in the inequality (23a). 

Some of the results of the preceding discussion become invalid for the 
singular case of @= 1/m. Physically this means that the force is applied at 
a point which would not move if the beam were to buckle laterally under the 
critical end moment mj (@#=7). This special case, which was briefly alluded 
to in the previous section, is treated in a general fashion in the Appendix. In 
its application to the illustrative example being discussed here, it is analyzed 
fully in what follows; a comparison with the Appendix shows that the general 
principles developed there are confirmed for the case under consideration. 

It can be verified easily that F(a), which is defined in Eq. (29), has no root 
in the range 0< a < 77 for the singular case of = 1/m. F(7) is indetermi- 
nate, but the customary limit procedure leads to a value of 1/z; hence, by 
Eq. (31), the boundary of the stable range is reached when w assumes the 
finite value governed by 


a, =7 (35) 
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while, by Eqs. (32), (33), and (34) and similar limit procedures 
Ky (1) = 
Kal?) = /- 
K =(/+ 7%) -! 
Similarly, the deflection mode approaches the limiting functions 
3 


(37) 


3 
2 


For values of w>wp, the restraint moment m retains its critical value 


mj, or a =m throughout. However, the solution (37) of Eq. (26) is no longer 
unique; it can in general be expressed by 


(38) 


in which c is an as yet arbitrary multiplier; the factor associated with c has 


been added for convenience, and (uy, By) is the normalized buckling mode given 
by 


U(Z)= cL) 4 sin =-2f,(2) (39) 


The value of c is now determined from the compatibility condition (27) 


which, in view of Eqs. (37), (38), and (39) and of the definitions of a and w , 
becomes 


=/-0 7 Yoru?) (40) 


Let a factor T $ 1 be defined by 


(41) 


where Wp is given in Eq. (35) and represents (in review) the value of w as 
a =7 is reached initially. Then c is related to T by the equation of the unit 
circle 
Cc 34 / (42) 

It is noted, and indeed expected from the discussion in the Appendix, that 
for any value of T< 1 there are two possible values of c; in particular, for 
T =0 (i.e., as the force A goes to infinity), c = +1. In other words, two dis- 
tinct equilibrium configurations are now possible, which are not adjacent to 


one another; thus, a “snap-through” (Durchschlag) phenomenon is to be ex- 
pected. 
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The reduction factors x are obtained by considering Eqs. (37), (38), (39), 
(41), and (42). This leads to the following set of relationships: 


Kg 2/772) 


(43) 


K 


I| 


-4 


In other words, the limiting configurations represent either predominant 
bending with little twisting, or else pure twist without any bending. Only one 
value of K appears, however, since the transition from one configuration to 
the others constitutes a rotation about the point of application of the force. 
Moreover, this value of x remains constant once the condition @ = 7, or 

W = Wo, is reached. For w > Wy (7 < 1), the force-displacement relationship 
is therefore linear, but the apparent stiffness of the structure is almost 
double that predicted by the linear theory. 

In concluding this section, it may be noted that, for sufficiently large values 
of the applied force, the response of the structure is a discontinuous function 
of the eccentricity e of the force. In fact, for “small” eccentricities 
(¢ < 1/7), the configuration approaches one of mostly bending and little 
twisting, the ratio being nearly independent of ¢. More surprisingly, per- 
haps, the response for “large” eccentricities (¢ >1/7) approaches one of 
pure twist without any bending, the reduction factors x, and kg being zero 
and unity irrespective of the value of ¢. For the singular case (¢ = 1/7), the 
two possible configurations discussed above represent limiting cases as the 
value of ¢ approaches 1/m from below or from above, respectively. 

This type of discontinuity seems to be one of the salient features which 
distinguishes the present theory from the conventional linear approach. 
Other examples have been investigated and have been found to lead to similar 
results. For example, if the same beam is subjected to two forces applied at 
equal distances from the ends and at equal (nonvanishing) eccentricities, but 
pulling in opposite directions, then the asymptotic response is found to be a 
discontinuous function of the ratio of the magnitudes of the forces. In fact, 
the larger of the two forces dominates the behavior entirely; singularity, that 
is, snap-through, occurs when they are equal. 


Experimental Results 


To obtain an experimental check on the results of the preceding section, a 
test arrangement similar to the one described(5) was employed. The beam 
specimen used was a strap, one inch high and 1/16 inch thick, which was 
made of heat-treated steel with a yield point of about 180,000 psi. Spanning 
a distance of.20 in., it was elastically restrained at both ends in the vertical 
plane, with the degree of end fixity y (see Eq. (27) ) computed at 0.74. How- 
ever, by loosening the clamps this end restraint could be removed entirely; 
in this fashion, the beam was made statically determinate to provide a check 
for the linear theory. The values for the elastic constants E and yu (Poisson’s 
Ratio) were known to be 30,000,000 psi and 0.3, respectively; this establishes 
the relationship A = 1.48w, in which A is the applied force measured in pounds, 
and w is defined in Eq. (30). 
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A relatively rigid vertical bar was attached to the beam at midspan; this 
bar was then subjected to a horizontal lateral force of increasing magnitude 
and varying eccentricity. The results corresponding to eccentricities e of 4 
in., 8 in., and 12 in. are discussed in what follows; these eccentricities are 
associated with values of ¢ equaling, very nearly, 1/27, 1/7, and 3/27, re- 
spectively, where the second value represents the singular case. The lateral 
deflections \u and rotations AB, as well as the total displacement A(u+ef) of 
the applied force, were measured in the usual manner by means of scales and 
mirrors. 

A comparison between the predicted and measured results is given in Figs. 
2, 3, and 4, which correspond to the three types of eccentricites (small, singu- 
lar, and large) mentioned above. In each case, the curves give the computed 
values of the reduction ratios on the basis of Eqs. (32),(33), and (34) as func- 
tions of the dimensionless loading parameter w (defined in Eq. (30) ). For 
the singular case of e = 8'', or ¢ = 1/z, the boundary of the stable domain is 
reached for wo = 4.06, as given in Eq. (35), with the associated reduction co- 
efficients computed in accordance with Eq. (36). For values of w in excess of 
Wo, these coefficients are governed by Eq. (43). 

The experimental points shown in the figures represent average values 
based on several test sequences. It is seen that reasonable agreement was 
obtained for small and for large eccentricities. Such quantitative discrep- 
ancies as do occur seem to be due to the effects of initial imperfections and 
of prestressing moments; for large values of w, other non-linear factors 
which were ignored in the present analysis cause further discrepancies. In 
fact, since both the measured deflections and rotations were somewhat in ex- 
cess of the limitations laid down in the Introduction, a more accurate analysis, 
which involves, for example, powers of u and 8 above the second, yields more 
compatible results. Nevertheless, there is excellent qualitative agreement 
between the theory and experiment; in particular, the discontinuous character 
of the asymptotic solution is fully corroborated. 

For the singular case (e = 8'', ¢ = 1/7) the agreement is confined to values 
of w below w,. For larger values of w, two possible solutions (Durchschlag) 
did indeed occur, but the quantitative agreement between the theory and the 
experimental values is poor. A possible explanation for this can be found in 
the effect of initial imperfections. These imperfections, such as initial cur- 
vatures, etc., become increasingly important as the boundary of the stable 
domain is approached, i.e., as the value of @ approaches 7. In fact, if Eq. (26) 
is properly modified to include the effect of an initial deviation from perfect 
shape, it can readily be verified that, in the general case, the solution “blows 
up” as the support moments approach their critical values. In effect, the 
orthogonality condition (A7) (see Appendix) is violated, which means that 
neutral equilibrium is approached only as the load goes toward infinity, in- 
stead of the finite value of w, predicted by the idealized theory. This may 
explain why the double-valued equilibrium configuration was “delayed” be- 
yond its theoretical value. 


CONCLUSION 


It may be stated that the theory presented herein constitutes a third ap- 
proximation to the problem of a beam subjected to a combination of bending 
and torsion. In the first approximation, the strains are assumed to be linear 
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functions of the displacements, while the equations of equilibrium refer to the 
undistorted configuration. This leads to an entirely linear theory, to which 
the principle of superposition applies. Such a theory, which is one commonly 
used in ordinary design, is valid in case the lateral stiffness of the beam is 
comparable in magnitude to its major stiffness. 

When this condition is violated, the need for a second approximation 
arises. In fact, the presence of large lateral deflections and rotations makes 
it imperative that the equations of equilibrium be referred to the distorted 
geometry. The resulting theory is still linear with respect to the lateral 
loads and twisting couples; however, the effects of the vertical loads can no 
longer be added linearly. 

The present, or third, approximation further abandons the assumption of 
linear strain-displacement relations. This is a natural result of the presence 
of large lateral displacements, in consequence of which the second approxima- 
tion appears to be illogical. This is not quite the case, however. Actually, at 
least so far as the lateral displacements and rotations are concerned, the re- 
sults of the last two approximations are identical for beams which are stat- 
ically determinate relative to moments in their major plane of stiffness. 

This may provide an explanation as to why this point has been ignored hereto- 
fore. In any event, it appears that, at least in the case of statically redundant 
structures, the introduction of nonlinear strain-displacement relations leads 
to substantial modifications in their predicted behavior. 
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APPENDIX 


_ In what follows it will be shown that, in general, the function U(M;u,8) 
= U(A,) (or, equivalently, W(A;) ) has an interior minimum, which is ap- 
proached asymptotically as the load parameter A goes beyond bounds. In 
addition, the singular case, in which such an interior minimum may not exist, 
is explored further. 

To this end, consider again the region of stability, which will be assumed 
to be bounded, as before. In the A,-space this region R is enclosed by the 
boundary B, on which the homogeneous equations (1) and (2) admit a 
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nontrivial solution (uy, By); on B the equilibrium is therefore “neutral,” with 
(uy, B;) representing the buckling mode. Moreover, in the general case, the 
inhomogeneous equations (1) and (2) admit no solution on B, while U grows 
beyond bounds as the stress point P (whose coordinates are the redundant 
parameters AY) approaches B. If now U can be shown to be convex, the 
existence of a unique interior minimum can be postulated. 

To show this, it is recalled, in line with footnote 7, that if u and # are dif- 
ferentiable functions of ),, then, by the inequality (18) or by an independent 
calculation, 


(Al) 


Consider now a plane which is tangent to U at some point (A,); its equation 
is given by 


er 


or, in view of Eq. (Al), by 
It follows therefore as a consequence of the inequality (17) that 
(A3) 

In other words, the surface U does not cross any of its tangent planes; it is 
therefore convex, and the existence of a unique interior minimum is proved, 
provided all the other previous assumptions are found to hold also. 

The exceptional (singular) case occurs when the inhomogeneous equations 


(1) and (2) do admit a finite solution. For the sake of simplicity, let Eqs. (1) 
and (2) be rewritten in the symbolic form 


Lu (4; 4 6) = 9 La(/g wu (inR) (A4) 
and let, in accordance with the previous discussion, 

Then the system of equations 

Luly, 4, B)= “s(4;% (on B) (A6) 


has a finite solution on that portion B' of B on which the “orthogonality con- 
dition” 


(on B') (A7) 


holds. Moreover, since U is finite on B', the existence of an interior mini- 
mum can no longer be postulated. 

For the purpose of studying singular behavior, let it be assumed that an 
interior minimum does not exist; this means in turn that U has an exterior 
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minimum, which lies on B'. In that event, as the load factor A increases, the 
stress point P reaches some point on B', from which it proceeds along B' 
toward the minimum. If B' is represented by a single point, as in the case of 
the illustrative example, the stress point stays fixed after reaching B' and 
the structure exhibits the quasi-linear behavior discussed earlier in the 
paper. 

After P has reached the boundary B', the most general solution of Eq. (A6) 
subject to the condition (A7) can be expressed in the form 


Cu, B) = (Ue, Bo) Br) (A8) 


in which (up, Bo) represents some solution of Eqs. (A6) and (uy,A4) is 
governed, subject to an arbitrary multiplicative constant, by Eqs. (A5). It is 
convenient to identify these functions further by the orthonormality conditions 


2 > 
678. 42 (A9) 


[(G? 46°67) =) 


in which b is an arbitrary number having the dimension of a length. Except 
for the sign of (uy, A), conditions (A9) determine all the functions uniquely 
for given value of b. 

Let the point at which the stress path “first” reaches the boundary B' be 
designated by P', and let P' be associated with a load parameter A‘ anda 
parameter c=c' (see (A8) ). Then the value of c' can be determined by a 
limit process as follows: 

While the stress point P is in the interior of the region of stability, the 
solution of Eq. (A4) can be represented by Eq. (A8), subject to the ortho- 
normality conditions (A9), if (uy, 84) satisfies the equations 


in which yp is the eigenvalue of the system (A10). This eigenvalue satisfies 
the relationship 


which is obtained by multiplying both sides of Eq. (A10) by uy and fj, re- 
spectively, by integrating, and by considering the second equation of Eq. (A9). 
Similarly, the condition 


(A11b) 
UolulM; Uy = 

( Bhp (4: AJ] fe) +8, Lg af) Meo 

is obtained, in which the second integral follows from the first through inte- 

grations by parts. Finally, from Eqs. (A4), (A8), (A11), and the linearity of 


the differential operators, the parameter c is seen to be given by the equa- 
tion 


c= th) (A12) 


As the point P approaches the boundary B, the eigenvalue yu approaches 
zero and the factor c increases generally beyond bounds. However, on B' 
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the numerator in Eq. (A12) also vanishes, as postulated in Eq. (A7); hence 
the fraction becomes indeterminate, and the following relationship is estab- 
lished by means of the usual limit considerations: 


C [94+ Az (A13) 
In Eq. (A13) and in what follows, a dot (such as in #1) signifies differentiation 
along the path, that is, Or " 
The denominator in Eq. (A13) is obtained from Eq. (Alla) through differen- 
tiation. This leads to 


= 6.) + B, Lg (M; 4z 


(A14) 
u,Lu( ™; t le, (M; uj, dz 


+ + 6 (kms a2 


The last integral on the right side of (A14) follows from the definitions of L, 
and Lg and from the fact that M=m, by Eq. (6a) 

The first integral in (A14) vanishes; this can be seen by substituting Eq. 
(A10) and by differentiating the second equation of Eq. (A9) with respect to Ag: 
The second integral vanishes similarly, since it can be converted into the 
first integral (and hence zero) through four integrations by parts and in view 
of the boundary conditions. Two more integrations by parts finally convert 
the third integral into 


J ms ( - dz (a15) 
The establishment of the numerator in Eq. (A13) proceeds along analogous 


lines. In fact, since (uy, By) satisfies Eq. (A4) at P' (for M=My}), it follows 
that the relationship 


tB,)dz = Le, Bo) + Ue, dz 


holds. The second equality in the preceding equation is the result of a num- 
ber of integrations by parts. However, if Eq. (A10) is differentiated with 
respect to Ag, and if uw is set equal to zero, it follows that 


La Bi) = Lal My + Ms Br)” 
La (4; =Le(M; B) + = 


on B'. When these relations are substituted in the last integral and the first 
equation of Eq. (A9) is considered, then, after some more integrations by 
parts, the relationship 
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J =- J Ms ( B, - 2K (A16) 


is established. Hence, by Eq. (A13), c' is expressed as the quotient of the 
two integrals appearing on the right side of Eqs. (A15) and (A16). 

Offhand, this quotient seems to depend on the direction of the approach 
path. That this is not the case can be seen from the following considerations: 

In the linear space in which the stress point P is defined, the moment 
component m, can in general be expressed in the form m, = % My, in which 
the set of 7, represents the direction cosines of the approach path relative to 
the coordinate axes. Similarly, the integrals on the right side of Eqs. (A15) 
and (A16) are representable as linear combinations involving these direction 
cosines. 

Now both denominator and numerator components vanish for all directions 
which are tangent to the boundary surface B'. Indeed, since 4 vanishes iden- 
tically throughout B, =0 for all directions other than the one normal to B. 
Similarly, the numerator is different from zero only in the normal direction, 
which can be verified by considering that the orthogonality condition (A7) is 
satisfied identically on B'. Hence the direction of the approach path is im- 
material, and the limiting coefficient c' can finally be expressed by 


(Ug" fy + uy "Bo 2K£,' dz (A17) 
2 frm - dz 


in which m, represents the moment component normal to B!. 

After the bounding surface of instability B' has been reached (for A= 3’), 
the stress point P travels on B' as the load parameter A is increased 
further. The response A(u,f§) continues to satisfy Eq. (A8); when this is 


substituted in the compatibility conditions (11), the following system of equa- 
tions is obtained: 


Cc! 


J Me kBi2) dz (A18t) 


The subscript n denotes the direction of the outer normal to B', while the 
subscript t designates the totality of all directions which are tangent to B' 
and for which the condition (A7) is identically satisfied. In the set of equa- 
tions (A18t) only one integral appears on the right side; two more integrals, 
similar in form to the second and third integrals in (A18n), can be shown to 
vanish. 

In general, P moves on B' in such a way as to minimize further the value 
of U+Up/ A2. With \ going to infinity, a limiting point Py is approached, 
which corresponds to the smallest value of U on B'. This point is found 
by setting the integrals appearing on the right side of the system of equations 
(A18t) equal to zero; a limiting value of c is then obtained by substituting the 
functions (u, 8) so determined in Eq. (A18n). (Note that there are two solu- 
tions.) 
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In special cases the stress point P may remain at P'. This happens when 
either P' is the only point on B' (such as in the illustrative example), or else 
when P' happens to be the point associated with the lowest value of U on B'; 
this latter case occurs when the vector P*P' is normal to the instability sur- 
face. In either case, the final state of stress designated by P' is generally 
reached when the load parameter is finite (A = A‘). For A > A’, the response 
functions A(uo, Bo ) increase linearly with A similarly to a statically determi- 
nate structure. “However, as before, two solutions for c are possible accord- 
ing to Eq. (Al18n). 10 In other words, ‘two distinct equilibrium configurations 
are possible such that one may be reached from the other only through a 
snap-through (Durchschlag) process. As before, a limiting value of c is ap- 
proached as the parameter A grows beyond bounds. 


10. In the (X,cA) plane, Eq. (A18n) defines a hyperbola. This becomes ap- 
parent from the fact that, if — A two, the right side can be con- 
verted into the form -\2W + cA2B + ¢2A2 fi, in which the dot represents 
differentiation in the n-direction. It is clear that 7 is negative, sincey 
is positive inside the region of stability and negative outside of it. 
Similarly, W is seen to be negative; indeed if it were positive then W 
would exhibit an interior minimum, in which case P would not reach (or 
stay) on the boundary. And finally, if W were zero, P' would be reached 
only in the limit. The quadratic form has therefore a negative discrimi- 
nant; hence it designates a hyperbola. 
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Discussion of 
“MINIMUM-WEIGHT DESIGN OF A PORTAL FRAME” 


by W. Prager 
(Proc. Paper 1073) 


Author’s Closure 


W. PRAGER,! M. ASCE.—The purpose of the paper was to investigate the 
influence of the usual linearization of the expression for the structural 
weight. For this purpose, it seemed sufficient to consider a single condition 
of loading. Minimum-weight design for multiple conditions of loading has 
been treated by Foulkes.* 


1. Prof. of Applied Mechanics, Brown University, Providence, R. I. 
J. Foulkes, “Linear Programming and Structural Design,” Proceedings, 


2nd Symposium on Linear Programming, National Bureau of Standards, 
vol. I, pp. 177-184 (1955). 
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Discussion of 
“WIND INDUCED VIBRATION OF CYLINDRICAL STRUCTURES” 


by Joseph Penzien 
(Proc. Paper 1141) 


ANATOL ROSHKO. !—The phenomenon of wind-induced vibrations of elas- 
tic structures is a very interesting problem in mechanics and one of great 
concern to the engineer. A basic understanding of the problem is still lack- 
ing. The main contribution of the present paper is to point out that strong 
vibrations of a cylinder may occur at the natural frequency even when this 
does not coincide with the vortex-shedding frequency. The author describes 
these as “self-induced” vibrations, as distinguished from those induced by 
vortex shedding. 

From the data given in the paper one can compute the Reynolds number for 
each velocity. It is seen then that the “self-induced” vibrations occurred at 
Reynolds numbers somewhat higher than 10° but lower than 106; that is, in 
the range of critical Reynolds number. In this critical range one ordinarily 
observes intermittent changes in the flow, between a configuration in which 
separation is upstream and one in which it is downstream of the 90° points of 
the cylinder. Since these movements of separation point are not necessarily 
symmetric with respect to flow center line, it may be expected that intermit- 
tent changes of lift will occur under these conditions. These randomly occur- 
ring changes of lift could produce the “bursts” of vibration at natural fre- 
quency, which the author describes as self-induced. Of course, the 
terminology may be quite appropriate, but it is important to know whether this 
phenomenon is related to the flow intermittency or whether it occurs also in 
other ranges of Reynolds number. 

The author’s presentation of the data would have been considerably more 
useful if the data had been reduced to coefficient form. As it is, one cannot 
see how the observed amplitude depends on the Reynolds number, the flow 
dynamic pressure and the dynamic characteristics of the elastic system. The 
reader may compute for himself the values of Reynolds’ number, from the 
data presented, but there is no way of calculating the actual aerodynamic 
forces or force coefficients, for the dynamic characteristics of the elastic 
system are not given, and were not taken into account in the data presented. 

The main criticism of the paper is that it is hardly more than a presenta- 
tion of the raw data, so that one can only speculate about the observed effects, 
and about their appearance under other conditions. 


1. Guggenheim Aeronautical Laboratory, California Institute of Technology, 
Pasadena 4, Calif. 
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Discussion of 
“BENDING OF A RECTANGULAR PLATE WITH ONE FREE EDGE” 


by W. Nachbar 
(Proc. Paper 1196) 


i. SILVERMAN, | A.M. ASCE.—With the drudgery of calculation removed 
by the use of a high-speed digital computer, the author has been able to pre- 
sent in numerical and graphical form the results of the analysis of a plate 
problem, the solution of which follows the classical methods of plate theory. 
The presence of a completely free boundary always makes the solution of a 
plate problem quite difficult, and it is for this reason that few attempts have 
been made to solve approximately the type of problem described in this paper. 
To construct approximate functions for w which satisfy the free boundary con- 
dition, as well as the supported boundaries, presents as difficult a problem as 
solving the equations and boundary conditions exactly. The difficulty lies in 
the moment and shear condition at the free edge which involves the value of 
Ww, W', WwW", and W''' at the free edge. This boundary condition can, of course, 
be simplified by assuming that the moment condition is Wh = 0 and the shear 
condition is W, " = 0, equivalent to assuming that u = 2 and uw = 0 
simultaneously. 

The Kirchhoff treatment of the plate problem furnishes conditions at the 
free edge which may be used in lieu of equations 20 (a,b). It can be shown for 
the problem treated by the author that the following relationships must hold: 


(N is used instead of the author’s D) 


J ( ) Swaxdy = 0 (A) 


onl Su Fu | (B) 
N | Qw x =0 
J 


1. Engr., Bureau of Reclamation, U. S. Dept. of the Interior, Denver, Colo. 
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(E) 


(F) 


(G) 


If w is taken as a product of two functions X,Y where Xo is chosen so that 
conditions (D) and (E) hold, then (A) will be satisfied if Y is a solution of the 
following differential equation 


a 
+ yre + q = fe) (H) 
The solution of (H) is given by 


Y =A ¥ily) +B Ve(y) + C ¥aly) + D valy) + f(y) (I) 


and the conditions (B), (C), (F), and (G) will be satisfied by appropriate choice 


4 
of the constants A, B, C, and D of (I). = 
The exaction solution of the problem satisfies equations (2) and 20(a,b) in Ec \f a 
lieu of equations (A), (B), (C) above. Ro 
The above sketched-out method has been used by the writer to solve plate Ene 


problems with free edges. Only slide rule accuracy is needed in most cases. 
For the plate shown in Figure 1, where the notation differs as shown from the 
author, under uniform load q, (K = 0) 


4 
(J) 
4 


Wily) = shay sin By 3 Vo(y) 


¥a(y) = shay cos By ; Waly) = chay cos py 0.026, 560 


a a 


The constants A, B, C, D of equation (I) are determined from the following 
simultaneous equations: 


ag¥''"(o) + va, = 0 (K) 


asY''*(0) + (2 azY'(o) =0 


(L) 


Y(h) =O (M) 


¥''(h) = 0 (N) 


The constants ag and ag are defined by 
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f(y) = an y where ¥ is the fluid unit weight so that ”h = qo. 
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Application to Author’s Case of 2 = 1; K = 0 (Uniform Load); v = 0.3 


Wi(h) = 0.6161 ; Wo(h) = 0.6184 ; Wo(h) = 12.5406; Wa(h) = 11.5037 


Equations (K) - (N) reduce to 


D = 0.4286K - 0.04552A 
C = 0.03151B 
0.0541A + 0.9825B + 5.964K = 0 
68.7A + 210B + 757.65K = 0 
A=+9.09KK ; B= -6.56K ; C = =-0.207%K ; D = = 0.0134K 


The deflection w(o0,o) which corresponds to the author’s w(o,b) is given by 
w(o,o) = 5(K + D) 


w(o,o) = 5 (0.9866)K = 4.933 4 ) = * = 0.140 


where ¢ = span and t = thickness of slab. 


(0,0) =-N x**(o)x(0) + 


(ox - )( 


using A = 9.09K; 0.3 


NK 
M,(0,0) = 10.82 & 


2 
(0,0) = 0-113 af 


aM,y(a,o) = = 2N(1 - v) 
2N(0.7) (- - ) (BB + Ca) 


0.854(6.56)4 


2 


0 .854(6.56) (4) = = 0.0583 
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It is seen that the proposed approximate method gives results which are coin- a 
cident with those of the author to at least two significant figures and which 43 { 
may be obtained by use of the slide rule. To indicate that such results are | 
not limited to the square plate, the table below compares coefficients for i 


w(0,0), M,(0,0), R(a/2,o) for various ( >) ratios, both for hydrostatic and 


uniform loads. The figures in parentheses are values taken from author’s 
tables and curves. 


TABLE 

Uniform Load (K = 0) al 
Coeff Coeff Coeff 
w(o,0) M (0,0) R(®/2,0) 


0.080 0.064 - 0.064 
(0.077) (0.060) (0.064) 
0.140 0.113 - 0.058 
(0.140) (0.112) (0.060) 
0.164 0.131 - 0.049 
(0.165) (0.132) (0.049) 


Hydrostatic Load K =1 


0.0248 0.0198 - 0.0177 
(0.0251) (0.0197) (= 0.0174) 
0.0397 0.0317 - 0.0087 
(0.0402) (0.0325) (- 0.0086) 
0.0328 0.0261 + 0.0009 
(0.032) (0.0256) (0) 


An iteration process is obviously possible since with 
the value of Y given by (I), a new value of X, may be 
obtained. 


Maximum departures from the exact solution occur when derivatives 
higher than the second are required, such as in the values for shear and edge 
reaction forces. For example, under uniform load the coefficients for 


vx(2,b) (author’s coordinates) are (0.16, 0.29, 0.33) to be compared with the 
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author’s values of (0.16, 0.25, 0.29) for ratios 0.5, 1.0, and 2.0. 

While comparable approximate results for uniformly loaded plates may be 
obtained by using a finite number of terms in the Fourier expansion, the pro- 
posed approximate method has been found useful in obtaining results for 
partial loadings, especially for free boundaries combined with conditions of 


fixity. 
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Proceedings of the American Society of Civil Engineers 


INELASTIC ANALYSIS OF ECCENTRICALLY-LOADED COLUMNS 


Marlyn E. Clark* A.M. ASCE, Omar M. Sidebottom* and 
Robert W. Shreeves** J.M. ASCE 
(Proc. Paper 1418) 


ABSTRACT 


This paper presents atheory by which the collapse load of an eccentrically- 
loaded column can be rationally determined. Experimental verification of the 
theory was obtained from tests made on aluminum alloy rectangular- and 
T-section columns. A study of the effect of time on the magnitude of the col- 
lapse load was also made using SAE 1020 steel columns. 


I. INTRODUCTION 


1. Preliminary Statement 


In the design of slender members subjected to compression loads, it is 
important to realize that the determination of the load corresponding to 
failure is usually a problem in unstable equilibrium. The failure is mani- 
fested in the form of a sudden buckling or collapse whenever the balance be- 
tween internal resisting moment and the externally applied moment is des- 
troyed. It should be noted that the determination of the failure load is not 
dependent on the critical stress level at which failure is initiated in the most- 
stressed fibers except for columns having small values of the slenderness 
ratio. Here the type of failure changes from one of instability to one of 
yielding or crushing. 

Stress values, in themselves, are of no concern whatsoever in the design 
of Euler columns—the size and shape of the cross-section of the members 


Note: Discussion open until March 1, 1958. Paper 1418 is part of the copyrighted 
Journal of the Engineering Mechanics Division of the American Society of Civil 
Engineers, Vol. 83, No. EM 4, October, 1957. 

* Associate Profs. of Theoretical and Applied Mechanics, University of 
Tllinois, Urbana, 

** Instructor in Theoretical and Applied Mechanics, University of Illinois, 

Urbana, 
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and the stiffness of the material are the only significant variables. For in- 
termediate values of the slenderness ratio, the collapse load is reached only 
after inelastic deformation has been initiated and its magnitude depends not 
only on the stiffness of the material but upon the shape of the stress-strain 
diagram for small inelastic strains. 

It can be said therefore, that an instability type of analysis must be used 
to calculate the failure load for all columns with slenderness ratios large 
enough to fail by buckling type of action. In such columns, the analysis must 
be capacle of relating the interval resisting moment with the externally ap- 
plied moment making use of the properties of the cross-sectional area and 
the stress-strain diagram of the material used. 

Several methods are available for designing centrally-loaded columns. 
When the material making up the column has a yield point, empirical equa- 
tions of many types have been used. For materials which strain-harden, the 
rational tangent or the double (reduced) modulus theories are available. 
When the column is eccentrically loaded, the maximum elastic load can be 
predicted by the secant formula, but the column will not collapse until inelas- 
tic action occurs. In some cases the collapse load may be twenty to thirty 
percent greater than the maximum elastic load, and this increase can be ob- 
tained by allowing only a small amount of inelastic deformation. 

The exact solution for determining the collapse load of an eccentrically- 
loaded column has not as yet been developed. By numerical procedures, the 


exact solution may be closely approximated, but from the standpoint of design, 


these solutions are extremely unwieldy. Various usable approximations to 


determine this load have been presented by Ketter,(1) Jordan,(2) Shanley,(3) 
and others. 


2. Purpose and Scope 


The primary purpose of this paper is to present a theory of inelastic 
column behavior which offers a relatively simple semi-graphical method for 
the determination of the collapse load. In the development of this purpose, 
analytical expressions for interaction curves* and moment-load curves** 
have been derived. The method utilized has been applied to rectangular- and 
T-section members; however, it can be applied to cross-sections of many 
other shapes. 

The theoretical value of the collapse load for an eccentrically-loaded 
column, as herein presented, is obtained as follows: first, moment-load in- 
teraction curves for various depths of yielding are constructed using the 
equations herein derived; second, the intersection of the theoretical moment- 
load curve with each of these interaction curves is located; finally, a smooth 
curve is faired in through these intersections and the maximum or collapse 
load selected graphically. 

The interaction curve for a given depth of yielding can be obtained by 
locating several points on the curve and then drawing a smooth curve through 


combinations of moment and load which, when applied to a column, result 
in a given depth of yielding or inelastic deformation. 


** A moment-load curve is a curve which, for a given loading system, relates 


changes in load with the corresponding changes in moment. 


* In this paper, an interaction curve is a curve which represents all possible 
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these points. To locate any one point, a strain distribution must be assumed, 
the stress distribution constructed using appropriate values from the stress- 
strain diagram, and the load P and the moment M calculated using the stress 
distribution and the equilibrium conditions. Five depths of yielding (0, 1/8, 
1/4, 1/2, and 3/4) have been considered so as to provide several points on the 
theoretical moment-load curve for a given column. 

The intersection of the theoretical moment-load curve for a column and the 
interaction curve for a given depth of yielding has been obtained by equating 
the radius of curvature of the column at the critical section to the radius of 
curvature needed to produce a given depth of yielding. In order to obtain an 
expression for the radius of curvature of the column, it is necessary to as- 
sume a configuration of the deflected axis of the column. Two assumptions 
have been made—first, an assumption of an arc of a circle, and second, the 
usual assumption of a cosine curve. 

An experimental program was undertaken to determine the ability of the 
proposed theory to predict collapse loads. Test data corroborating the theory 
were obtained from rectangular columns made of annealed SAE 1020 steel 
and rectangular- and T-section columns made of 24ST aluminum alloy. 

A secondary purpose of this paper is to study the effect of time on the 
magnitude of the collapse load. Since previous investigations on mild steel 
beams by Clark, Corten, and Sidebottom(4) indicated that the yield strength of 
mild steel can be lowered if sufficient time is allowed for inelastic deforma- 
tion to be completed at a given load, the steel columns were tested under 
similar conditions in a constant-load testing machine.(5) Several columns 
were loaded in steps and each successive load was maintained for 6 minutes 
or longer. The results of tests under these conditions indicate that the yield 
strength was from 5 to 10 percent below that usually obtained in an ordinary 
tension test. This is of the same order of magnitude as the value found in the 
beam investigation.(4) Other mild steel columns were tested so as to com- 
plete the loading to failure in less than two minutes. These data indicate that 
the yield strength of the material in the column had been raised in the 
presence of rapid loading. 


tl. Theory 


3. Derivation of Moment and Load Interaction Equations 


The derivation of interaction equations for a given depth of yielding is 
based on the assumption that plane sections in the column remain plane and 
that the stress-strain diagram of the material may be represented by two 
straight lines as indicated by Fig. 1. Depth of yielding has been chosen as 
the invariant (instead of some other quantity such as maximum strain or 
maximum stress) because of the ease with which the derivations are made 
when this quantity is chosen. 

The interaction curve for any given depth of yielding is usually made up of 
one or two linear segments and one curved segment. The interaction curve is 
linear so long as yielding has not occurred on both sides of the member.(6) 
Therefore, the most convenient method of constructing an interaction curve 
for a given depth of yielding is to (1) locate the terminal points of the linear 
segments, (2) determine the coordinates of several points on the curved por- 
tion, and (3) “fair in” a smooth curve. 

To illustrate the method of computing a point on an interaction curve, 
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consider the column of rectangular cross-section shown in Fig. 2. The 
column has a width b and a depthh. For the strain distribution, the depth of 
yielding on the compressive side of the column is ay and the depth of yielding 
on the tension side of the column is a2. 

The stress distribution shown in Fig. 2 was constructed from the strain 
distribution and the stress-strain diagram for the material as shown in Fig. 1. 
It will be noted that the stress-strain diagrams for tension and compression 
are represented by two straight lines. In the elastic portion, the slope of the 
straight line is equal to the modulus of elasticity, E. In the inelastic portion, 
the slope is for compression and for tension where and may 
be considered the strain-hardening factors for the material. 

The load P and the moment M on the column may now be determined by 
making use of the stress distribution K-L-C-A-H-F-J. The work involved in 
using the equations for P and M will be reduced if they are derived so that P 
and M are functions only of aj, the depth of yielding on the compression side, 
ag, the depth of yielding on the tension side, and a, the elastic depth or core. 

In calculating the load P, it should be noted that for elastic conditions 
(such as stress distribution K-L-D-G-J), the load P would be equal to the 
product of the cross-sectional area A and the stress at the centroid 7,''' 
The desired load P is equal to this load for elastic conditions minus the load 
corresponding to the stress distribution A-C-D plus the load corresponding 
to the stress distribution H-F-G. From similar triangles, it can be seen that 


where the yield stress o, is the average of 0,1 and 02. The stress at the 
centroid then becomes 


«esa (c ) (1) 
& 1 


To obtain the loads corresponding to the stress distributions A-C-D and 
H-F-G, it is necessary to obtain the stresses represented by C-D and F-G. 
Stress C-D becomes 


20 


(l-@ )(o'-o )s(l-a) (2) 
a 


and an identical expression (with suitable subscripts) gives stress F-G. Us- 
ing the relations given in Eq. 1 and 2, an expression for the load P may be 
written as 


ho. be be" 


(3) 


Ps 29, be - 


In calculating the bending moment M corresponding to stress distribution 
K-L-C-A-H-F-J, it is convenient to subtract from a fictitious elastic moment 
M] (equal to oyI/c from the stress distribution K-L-D-G-J) the moments of 
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the loads resulting from stress distributions A-C-D and H-G-F. Using simi- 
lar triangles, the elastic moment can be written as 20,.I/a. The bending mo- 
ment M can now be written as shown below. 


20,1 ba? o ba® 


(2 a) (e - - (ec - 2) 


M= 


= - 


It should be noted that Eq. 3 and 4 are valid if either aj or ag is zero. 

Attention will now be turned to the derivation of the load and moment rela- 
tions for a T-section member. The cross-sectional dimensions of this mem- 
ber are shown in Fig. 3. The calculations are more complex than those for 
the rectangular section because of the change in width of the section. The ex- 
pressions for the load P and the moment M will be different depending upon 
whether the depth of yielding on the compression side, aj, is less than or 
greater than the thickness of the flange, tg. If aj is less than or equal to t2 
the expressions for the load and moment are similar to Eq. 3 and 4. Taking 
into consideration the difference in width of the section for the stress distri- 
bution A-C-D and H-F-G, the load and moment are given as follows: 


Peo, A- +(1-a) == 
20 I o ba® 
Me (1-4) -=) 6) 


a a 2 


If the depth of yielding on the compression side of the T-section member 
extends through the flange, i.e., if aj is greater than tg, the calculations of 
the load and moment for the stress distribution A-C-D is complicated by the 
abrupt change in width of the section. The computation of the load can be 
simplified by calculating the load corresponding to the area ajb and then sub- 
tracting from this the load corresponding to area (aj - tg)(b - t). From 
geometry, the magnitude of the stress at the inside of the flange minus the 
yield stress o,, is given by the relation 20¢ (aj - tg)/a. Using this relation 
along with Eq. 1 and 2, the load expression becomes 


ey a a a 


20_A (b-t)(a -t ta? 
Peo, - -0, —— 4(1-0 ) (7) 


The determination of the moment corresponding to stress distribution 
A-C-D is simplified if it is realized that the moment of the stress distribu- 
tion A-C-D is equal to the product of (1 - 1) and stress distribution A-B-D. 
But this calculation is complicated by the change in width of section over the 
depth aj. This difficulty can be circumvented by noting that the magnitude of 
the fictitious elastic moment My depends only on the stress gradient since a 
uniform stress distribution superimposed on a bending stress distribution 
does not affect the magnitude of the resisting moment. Therefore, the mo- 
ment resulting from the stress distribution A-B-D is equal to Mj minus the 
moment of the stress distribution A-E-G which does not have a change in 
width. Thus, 
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+ a_)*t 


ata 
M= M,-(1 a) M, = (la) 


Again it is noteworthy that aj and a2 can assume the value of zero without 
affecting the validity of Eq. 5, 6, 7 and 8. 


4. Interaction Curves 


Interaction curves for a column of rectangular cross-section have been 
constructed by the use of Eq. 3 and 4 and are shown as the solid curves in 
Fig. 4. These curves are for a mild steel having equal yield stresses in 
compression and tension, i.e., 0g1 = Ge2 = Se, and for a strain-hardening 
factor @ equal to zero. It will be noted that for a depth of yielding not equal 
to zero, the interaction curve is made up of a curved portion and a straight 
line. For instance, consider the curve A-B-C for three-quarter depth yield- 
ing. From A to B, yielding occurs on both faces of the member. From B to 
C, yielding occurs over the same total area of cross-section, but the stress 
on the compression side remains constant at 0, while the stress on the ten- 
sion side of the member varies from 0, tension to Se compression. The 
proof of the linearity of line B-C may be found in a paper by Brush and 
Sidebottom. (6) 

For the T-section columns, the interaction curves can be obtained by mak- 
ing use of Eq. 5 through 8. Typical curves for the aluminum alloy columns 
are shown in Fig. 10 to 12. It can be noted that all interaction curves have 
been plotted in dimensionless form. The load and moment equations have 
been divided by the maximum elastic load, Pa, and the maximum elastic bend- 
ing moment, Me, respectively. In these expressions the compressive yield 
stress 0, should be used to determine Pe, and either Oe or de2 should be 
used to determine M, depending upon which yield stress is reached first when 
the member is subjected to pure bending. Accordingly, cj or c2 should be 
selected to correspond with the appropriate stress. 


5. Development of Moment-Load Curves 


A moment-load curve for a column is a plot of corresponding values of the 
load P and the moment M. The moment M can be seen to be directly depend- 
ent on the eccentric load and equal to PA, where A is the initial eccentricity 
e plus the deflection 5 of the column. It is the stated purpose, therefore, of 
this section to set forth a procedure whereby the moment-load curve can be 
determined analytically. Once the moment-load curve for a given column is 
known, it is a simple matter to determine the collapse load from it. The load 
corresponding to a point on the moment-load curve at which an increase in 
load is impossible can be defined as the collapse load, P,. This point can be 
found graphically by drawing a vertical line tangent to the moment-load curve. 
The abscissa of this point of tangency is the ratio of the collapse load to the 
maximum elastic load, Pe. The method herein presented makes use of the 
interaction curves for determining the moment-load curve. 

Instability and collapse of an eccentrically-loaded column will take place 
when the externally applied moment exceeds the internal resisting moment 
offered by the type of material and the shape of the cross-section of the 
column. Stability can then be assured if there is a balance between these two 
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moments. A moment-load curve for stable conditions can be derived on this 
basis as follows: The interaction curves shown in Fig. 11 through 15 give 
combinations of axial load and moment which may be applied to the member 
to produce a given depth of yielding. For a column with a given eccentricity 
and deflection pattern, only one point on each curve represents the internal 
resisting conditions which prevail in this loading. This one point, therefore, 
is a unique point on the moment-load curve for this depth of yielding; when a 
method is found to determine the intersection of the moment-load curve with 
each interaction curve, it will be possible to locate enough points so that a 
representative curve can be drawn. 

It should be noted that each interaction curve is a curve for a constant 
depth of yielding. For the non-linear portion of a given curve, yielding oc- 
curs on both sides of the member and the radius of curvature of the member 
remains constant. When the expression for the radius of curvature of the 
column, as determined from its assumed deflected axis, is equated to the 
radius of curvature expression, as determined from the interaction curve, it 
is equivalent to equating the applied moment to the resisting moment (since 
moment and curvature are proportional if the deflection is small compared 
with the length of the column). In order to determine the radius of curvature 
for a column, the configuration of the column must be known. Since the exact 
configuration is impractical if not impossible to obtain, one must be assumed. 
In this paper, two assumptions will be presented. The first assumption will 
be that of a segment of a circle—this will yield a conservative estimate of the 
collapse load since each section of the column is assumed to be stressed the 
same as the central section. To provide a means of comparison with work 
done in this field by other investigators, the second assumption will be that 
of a cosine curve configuration. 

By assuming that the column deforms into a segment of a circle, the radius 
of curvature can be obtained from the expression (see Fig. 5) 


R? = (2)? + (R - 8)? 


For the cosine curve assumption, an expression involving the radius of curva- 
ture at the center of the member can be obtained if the following equation is 
differentiated twice with respect to x. 


mx 
= 60cos— 
4 L 


Assuming small deflections, the following simple relations between radii of 
curvature and deflection are obtained. 


2 
(circle) (cosine) (9) 
n 


The radius of curvature can also be found from the strain distribution in 
the column. For the strain distribution shown in Fig. 2, the strains at both 
edges of the elastic material are known and are equal to €,j and €¢2, the 

yield strain in compression and tension, respectively. Since plane sections 
are assumed to remain plane throughout the straining, a consideration of this 
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geometry of deformation will show that the reciprocal of the radius of curva- 
ture can be written in terms of these strains and the distance, a, between 
them. When the elastic core distance, a, is set equal to kh (where k is a 


proportionally constant denoting the elastic core depth), the expression can 
be written 


(10) 


where €, is the average of the yield strains € .j and €e2. It should be noted 
that k remains constant for any given depth of yielding. By equating Eq. 9 
and 10, the deflection of the column can be obtained in terms of the geometry 
of deformation, the dimensions of the cross-section, and the properties of the 
material as follows: 


aL? 
5 = (circle) : (cosine) 
Lich n= kh 


(11) 


Eq. 11 is valid as long as yielding occurs on both faces of the member. 

By making use of Eq. 11, it is now possible to construct the moment-load 
curve for the column for each assumption. The moment M can be written in 
terms of the load P, the eccentricity e, and the deflection 6. Using the circu- 
lar segment representation, 


4kh 


In order to make use of the dimensionless interaction curves, the above ex- 
pression must be divided through by the maximum elastic moment Me (equal 
I 


~ , whichever is valid) while the right side of the expression 


c 
2 
must be multiplied and divided by the maximum compressive load P, (equal 


0,11 
to eA). The expression which results, if Me equals - , is 


1 


By dividing by bo , the above relation can be put in a form useful for enter- 
e 


ing the interaction curves. This result is 


M/M, 


P/P, 


where p 2 = I/A. In Fig. 4, the line OD of slope @ was obtained for 
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three-quarter depth yielding by substituting k = 1/4 in Eq. 12. The intersec- 
tion of this line and the curve for three-quarters depth of yielding is point D, 
a point on the moment-load curve O-L-D. By assuming other values for k, 
other intersection points between the interaction curves and the moment-load 
curve are obtained. 

Since Eq. 12 is only valid when yielding occurs on both sides of the mem- 
ber, it is not valid for the straight line portion of the interaction curve. This 
difficulty is easily overcome since the deflection of the column at any point 
on the straight line portion of the interaction curve is proportional to the 
bending moment. The proof of this statement is given in Appendix A. Thus, 
the deflection corresponding to point L in Fig. 4 is 


M (13) 
8, 
L ty 


where 6 y is the deflection of the column corresponding to point U the end of 
the straight-line segment in Fig. 4 and is obtained from Eq. 11. My is the 
value of the moment at this point. My is the value of the moment at some 
lower point where the moment-load curve intersects the interaction curve. 
Since My./My appears in the equation for tan @ in the case when there is 
yielding on one side only, Eq. 12 must be solved by trial and error. For 
every value of My that is selected, there is a corresponding value for tang. 
When the value chosen for tan @ is equal to the right hand side of Eq. 12, the 
desired point on the moment-load curve is determined. The following equa- 
tions summarize the foregoing: 


(14) 


for M 
e 


ceo 
2 > 
fe) kh en My 


e2 


(16) 


2 


If the cosine curve configuration was to be utilized instead of the circle, Eq. 
14 through 16 would be identical with those given above except that the num- 


2 
ber 4 would be replaced by "> ‘ 
Il. Experiment 


6. Description of Materials and Members 


In the experimental program to determine the correlation between theory 
and experiment, a total of 4 columns made of 24ST aluminum alloy and 9 
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columns made of SAE 1020 steel were tested. The aluminum alloy was 
chosen to represent those materials which strain-harden at the onset of in- 
elastic deformation and which are nominally time-insensitive. The steel was 
known to possess time-dependent characteristics, and, therefore, it was 
chosen to study the effects of time (rate of loading and time at a given load) 
on the inelastic behavior of eccentrically-loaded columns. 

Two different types of cross-sections (rectangular- and T-section) were 
used for the aluminum alloy columns whereas only rectangular columns were 
made from the steel. In the as-received condition, both materials possessed 
a variation in properties across the diameter of the bar-stock that was used. 
A survey of this variation (made using plate tension specimens 1/4 x 1/2 in. 
in cross-section and 10 in. long) showed that the yield stress for both materi- 
als was approximately 10 percent greater at the periphery of the bar than at 
the center. 

Fig. 6 shows the three different types of column specimens used, their 
dimensions, and their locations in the original bar-stock. In Fig. 7 are shown 
three different sets of stress-strain data for the aluminum alloy material. 
Because of the variation in properties mentioned above, in each set of curves 
in Fig. 7 the properties reported were obtained from specimens cut from the 
bar stock near the place where the significant stress occurred in the column. 
For example, in the T-section column the critical tension fibers were located 
near the periphery of the bar, hence, the tension specimens were cut from 
this location. For the steel columns, the variation in properties was account- 
ed for by a greater number of compression tests, many specimens being cut 
directly from the undeformed ends of the columns. An overall average of 
30,000 psi, obtained from some 24 yield point values from both tension and 
compression tests, was used. 

In order to be utilized in the theoretical relations developed in Section II, 
it is necessary for the stress-strain diagram to be represented by two 
straight lines. No difficulty is presented for the SAE 1020 steel since the 
yield range is well represented by just a horizontal line. However, the dia- 
grams for the aluminum alloy shown in Fig. 7 present more of a problem. 
The question arises as to whether the straight line in the inelastic region 
should represent the stress-strain diagram accurately for large inelastic 
strains or for small inelastic strains. In general, the type of columns tested 
in this investigation reached their collapse load before excessive deformation 
took place and the small strain representation should be more apropos. Both 
representations were investigated and are discussed in Section IV. 

The accepted published value of the modulus of elasticity were used for 


both materials; 10,600,000 psi for the 24ST aluminum alloy and 30,000,000 psi 
for the SAE 1020 steel. 


7. Method of Testing 


Fig. 8 shows a sketch of the double knife-edge column fixtures and Fig. 9 
shows a typical set-up in the constant-load testing machine. Note that the 
two knife-edges are placed pernendicular to each other and point loading is 
achieved. The eccentricity was pre-set by moving the column seat to one 
side or the other of the center line of the knife-edge. All rectangular-section 
members were supported in the manner shown in Fig. 8; however, the T- 
section members were supported by only one knife-edge at each end (thus, 
line loading was achieved) in order to insure that the column would buckle 
about the desired axis. Since the effective length of a column must be 
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measured between the knife-edges, a distance of 2.50 in. (1.25 in. for each 
end of the column) was added to the length of each machined column to ac- 
count for the length of the knife-edge fixtures. The aluminum alloy columns 
were provided with increased dimensions of the ends so as to avoid local 
crushing of the material. 

The deflection at the center of each column was obtained by a dial mounted 
on the fixture shown in Fig. 9. It should be noted that the fixture was sup- 
ported by the blocks which contained the knife-edge seats. Provision was 
made for a deflectometer, capable of generating an electrical impulse pro- 
portional to the deflection, to be connected in series with the dial so that a 
continuous record of the deflection could be made on a Sanborn Recorder. 

The eccentricity for each column was adjusted on the fixtures and the mag- 
nitude was checked analytically using measured values of the maximum ten- 
sion and compression strains (obtained from SR-4 strain gages mounted on 
the column) and measured values of the deflection. The computed value cor- 
related well with the value set on the fixtures. The accuracy of SR-4 strain 
gages was checked in bending before the columns were tested under eccentric 
compression loads. 

Most of the columns were tested in a special constant-load testing machine 
which was designed and built for this investigation. The machine is described 
in a paper by Clark and Sidebottom.(5) The load was measured by a 20,000 Ib. 
load cell in which SR-4 strain gages were utilized for indicating load. The 
load was applied to the column by means of a hydraulic jack. The jack was 
connected to a ballast-type reducing valve through an oil-gas accumulator. 
The reducing valve maintained a constant pressure on the hydraulic system 
(and hence, a constant load on the column) by continually adjusting the pres- 
sure to compensate for any drop in pressure (and load) due to yielding in and 
deflection of the column. The reducing valve was connected to a high pres- 
sure source of nitrogen which it reduced to the correct value and introduced 
into the system when a drop-off in pressure was indicated. 

All the steel columns except Column SR7 were tested in the special con- 
stant-load machine. In several of the tests on steel columns, a step method 
of loading was adopted. In this method, when the stresses reached the inelas- 
tic range, each successively larger load was maintained on the column for a 
given period of time unless indications were that collapse was imminent. In 
this case the load was maintained until the column collapsed. Six-minute 
time intervals were selected for use in this series of tests. In one column 
test which started out with six-minute load intervals, the collapse load was 
maintained for approximately an hour before final collapse occurred. Three 
of the columns were in the inelastic range for a total time period ranging 
from 6 sec. to 82 sec. For these columns the load, deflection, maximum com- 
pression strain, and maximum tension strain were recorded automatically on 
a four-channel Sanborn Recorder. 

One of the steel columns (SR7) was tested on the 20,000 lb. range of a 
200,000 lb. Riehle screw-power testing machine in order to determine data 
beyond the collapse load. The screw-power testing machine was equipped 
with a load-holder, making it possible to maintain the load on the column as 
with the constant-load testing machine. 

The aluminum alloy columns AT1 and AT2 were tested in the Riehle test- 
ing machine, and columns AR1 and AR2 were tested in the constant-load test- 


ing machine. The schedule of loading was similar to that for the steel 
columns. 
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IV. Analysis of Results 


8. 24S-T Aluminum Alloy Columns 


The test results for the T-section columns are presented in Fig. 10, 11, 
and 12, and the curves for the rectangular-section columns are shown in 
Fig. 13. Pertinent information for all columns is presented in Table 1. 

The T-section columns had cross-sectional dimensions approximately 
equal to those given in Fig. 6. The thickness t in Columns AT1 and AT2 was 
0.24 inches and 0.25 inches, respectively, and the eccentricity of load for 
each column was 0.230 inches and 0.259 inches, respectively. The experi- 
mental moment-load curve for Column AT1 is shown in Fig. 10. The applied 
moment was computed at various stages of loading using the value of the de- 
flection as measured at the center of the column and the pre-set eccentricity. 
The theoretical circular segment and cosine curve approximations to the 
moment-load curve were obtained by making use of Eq. 7 and 8 (to determine 
the interaction curves) and the procedure outlined in Article 5. These curves 
are shown in Fig. 10 by the dotted and dashed lines, respectively. 

In constructing the theoretical curves shown in Fig. 10, it was necessary 


that each of the tension and compression stress-strain diagrams of the mate- "| 
rial be represented by two straight lines. Since the stress-strain diagram ES 3 
for 24S-T aluminum alloy is curved at the region of transition from elastic i= 


to inelastic behavior, it was impossible to find two straight lines which would 
represent the stress-strain diagram for all inelastic strains. Therefore, two 
idealized stress-strain diagrams were considered—one for the case in which 
small inelastic strains were closely approximated (this hereafter will be re- 
ferred to as Case I), and the other for the case in which large inelastic 
strains were closely approximated (hereafter Case II). Actually, Column AT1 
collapsed at measured strains of about 0.0065 in both compression and ten- 
sion; therefore, the best agreement should be seen for Case I (see Fig. 7). 

In order to gain some idea of how critical the representation of the stress- 
strain diagram might be, the data for Column AT1 were analyzed using both 
Case Iand II. Accordingly, the interaction and moment-load curves in Fig. 
10 have been computed making use of the stress-strain diagram as represent- 
ed by the dotted lines in Fig. 7, while for Fig. 11, the dashed lines in Fig. 7 
were used. It can be seen that the shapes of the curves are greatly altered by 
the two different representations. Regardless of this great difference, how- 
ever, the test data for Column AT1 lie between the two approximations of the 
moment-load curves in both Fig. 10 and 11. 

As indicated in Table 1, the test value of the collapse load for Column AT1 
was found to be 9,840 lb. The magnitude of the collapse load based on the 
circular segment approximation is 9,520 lb. and 9,740 lb. for Case I and II, 
respectively. Corresponding values for the cosine curve approximation are 
10,750 lb. and 10,450 lb. The circular segment approximation is based on the 
assumption that each section of the column is subjected to the same moment 
as the critical central section; since this is not the case, the theoretical 
deflection will be larger than the actual deflection and the theoretical collapse 
load will be smaller than the actual collapse load if the stress-strain diagram 
has been accurately represented. Fortunately, the representation of the 
stress-strain diagram is not too critical for determining collapse loads. The 
theoretical collapse load for Column AT1 (9,740 lb. by the circular segment 
approximation) is less than the actual value even if use is made of the 
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stress-strain diagram for Case II (i.e. for strains of about 2 per cent). This 
is a rather poor representation since the strains at collapse were less than 
1 per cent; however, this did cause a shift of the theoretical load toward the 
actual value. Similar conclusions could be drawn from the results of other 
trial representations of the data for the aluminum alloy columns although 
those data are not presented. 

The test data for the T-section Column AT2 are presented in Fig. 12. The 
theoretical curves were obtained using Case I. The results are similar to 
those found for Column AT1 in Fig. 11. It can be noted that for Columns AT1 
and AT2, the experimental moment-load curves show data beyond the col- 
lapse load. This was made possible by using the screw-powered testing ma- 
chine. For all other columns except SR7, the data cease at the collapse load 
Since at this point (under conditions of a constant load) catastrophic deflec- 
tion took place. 

The rectangular-section columns were tested in the special constant-load 
testing machine. The dimensions of the two columns were identical (see Fig. 
6) and the eccentricity for each column was 0.32h (0.259 in.). The data for 
both columns are plotted in Fig. 13. The theoretical curves were plotted us- 
ing Case I for the stress-strain representation. The test data indicate that 
the two columns were almost identical in behavior. Again, the experimental 
collapse load was found to lie between the two theoretical approximations 
(see Table 1). 


9. Mild Steel Columns 


The mild steel columns were tested primarily to study the effect of time 
on the heterogeneous type of inelastic deformation found in mild steel materi- 
al. Results of previous investigations(4) indicated that inelastic deformation 
continued for long periods of time at constant load and that the moment cor- 
responding to a plastic hinge was approximately 10 percent less than the cal- 
culated hinge moment (based on conventional stress-strain properties of the 
material) when the tests were carried out in a dead-load testing machine. In 
this same investigation, the yield point, as determined from tension members 
subjected to the same schedule of loading as were the beams, was found to be 
the same as that computed from the hinge moments in the beams. Thus, in 
the column investigation, it would be reasonable to assume the best compari- 
son of theory and experiment would be obtained if the stress-strain specimens 
were subjected to the same schedule of loading as the columns. This would 
be impractical for design purposes, however, because the stress-strain data 
as commonly reported are obtained from conventional testing procedures. In 
view of this, the tension and compression specimens were tested at a rate of 
loading such that the strain at the end of flat portion of the stress-strain dia- 
gram was reached in from three to five minutes. 


a) Sustained Loading Tests. All of the mild steel columns, except one, 
were tested in the specially designed constant-load testing machine. Since 
time at a given load and rate of loading were both believed to have a consider - 
able influence on the collapse load, six columns were loaded in a stair -step 
manner with each load being maintained for a period of time. The increase in 
each increment of load was from two to three percent of the yield point load 
and the load was maintained for six minutes unless the inelastic deformation 
continued to increase as it did at the collapse load. In Table 1 are listed the 
times required for each column to fail at the collapse load. In one case 
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inelastic deformation continued for almost an hour before final collapse. 

The results of the sustained load column tests are shown in Fig. 14 and 
listed in Table 1. It will be noted that the collapse loads for the columns sub- 
jected to sustained load all fall below the theoretical value predicted by the 
circular segment representation. The theoretical curves were based on a 
lower yield point stress of 30,000 psi as obtained from tension and compres- 
sion tests. Theoretically, the circular segment representation should be con- 
servative; however, the average collapse load was about 6 per cent below this 
value. This would indicate that the actual lower yield point stress for the ma- 
terial was at least 6 percent below the value obtained from the conventional 
tension and compression test. Thus, the time-dependent heterogeneous type 
of yielding which is characteristic of mild steel was influential in bringing 
about the same decrease in load-carrying capacity for columns subjected to 
sustained loads as for beams subjected to sustained loads. Furthermore, the 
collapse load data show considerably more variation than those for strain- 
hardening materials such as the aluminum alloy. 


b) Short-Time Loading Tests. Three columns were loaded to failure 
through the inelastic range in time intervals of 6, 20, and 82 sec. The results 
shown in Fig. 15 and Table 1 indicate that the collapse load has been appreci- 
ably increased. Since the greatest increase resulted for the column with the 
fastest rate of loading, these tests indicate the dependence of the collapse 
load on time. It is interesting to note that the variation in the yield point with 
the rate of strain is in itself sufficient to account for the variation in the mag- 
nitudes of collapse load. These three columns collapsed at a strain of ap- 
proximately 0.40 percent. The strain at the beginning of yielding was 0.10 
percent. Therefore, the rates of strain for Columns SR9, SR11, and SR15 
through the inelastic region were approximately 15 x 10-5, 50 x 1575, and 
3.7 x 10-° per sec., respectively. The strain rate used in the testing of ten- 
sion specimens was approximately 5 x 10-5 per sec. Thus, Columns SR9 and 
SR15 were tested at about the same strain rate as the tension specimens and 
the experimental collapse loads for these two columns fall between the two 
theoretical values, i.e., the theory is corroborated by the test data. Column 
SR11, which was loaded in 6 sec., had a strain rate that was about one order 
of magnitude greater than that used in the tension testing. It has been found 
in rate of loading tests on mild steel tension specimens ) that, in this range 
of strain rate, an increase of one order of magnitude will cause an increase 
in yield point of about 7 percent. If the yield point for Column SR11 had been 
increased 7 percent, the collapse load for this column would have fallen be- 
tween the theoretical values. More data are needed to strengthen this con- 
clusion, however. 

It should be noted in Fig. 15 that there is an increase in the moment at 
which the column collapses with increased speed of testing. For Columns 
SR9 and SR15, the collapse loads as well as the depth of yielding at collapse 
were found to confirm the theoretical analysis. However, the load on the 
column was found to remain constant for a considerable increase in moment 
while the theoretical curves show a decrease in moment as soon as the col- 
lapse load is reached. This difference can be attributed to the time-sensitive 
behavior of the mild steel. The heterogeneous yielding in the form of 
Lueder’s bands requires a certain length of time to be initiated and to be 
completed. Hence, any increase in the rate of straining will cause an in- 
crease in the collapse load and will also cause this collapse load to remain 
constant for greater deflections while the yielding is taking place. 
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. SUMMARY AND CONCLUSIONS bs 
10. Summary 
This investigation was conducted to obtain a more simplified procedure for > 
determining the collapse load of an eccentrically-loaded column. The appli- a 


cation of this procedure to columns made from time-sensitive and time- 
insensitive materials was also studied. The simplified procedure involved 
the use of constant depth of yielding interaction curves. Also involved was 
the construction of moment-load and load-deflection curves. These curves 
necessitate an approximation of the deflected axis of the column. Two ap- 
proximations were made—(1) the column was assumed to deflect into a seg- 
ment of a circle, and (2) it was assumed to deflect into the shape of a cosine 
curve. 

Rectangular- and T-section 24S-T aluminum alloy columns were tested to 
determine the validity of the approximations used in the theory since this 
material is known to be relatively time-insensitive. Nine rectangular-section 
SAE 1020 steel columns were used to investigate the effect of the time- 
sensitive behavior of this material on the collapse load. Six of these columns 
were tested in a constant-load testing machine where each inelastic load was 
maintained for six minutes unless inelastic deformation continued to increase 
as it did at the collapse load. The time required for collapse at these sus- 
tained loads was as much as one hour. Three steel columns were also loaded 
through the inelastic range in time intervals varying from 6 sec. to 82 sec. 


11. Conclusions 


1. Once the interaction curves for columns of a given material and cross- 
section are known, the theoretical moment-load and load-deflection curves 

for the column can be quickly constructed. The collapse load then can be ob- 
tained from either curve. 


2. The theoretical equations developed require that the stress-strain dia- 
gram be represented by two straight lines. For aluminum alloy material, the 
transition from the elastic to the inelastic region of the stress-strain dia- 
gram is somewhat uncertain. Therefore, the experimental collapse load data 
were analyzed using two different representations of the stress-strain dia- 
gram (see Fig. 7). Both of these representations were considered in the con- 
struction of the interaction curves and the theoretical moment-load curves. 

It was found that the representation of the stress-strain diagram greatly in- 
fluenced the shapes of the theoretical curves but that the effect on the col- 
lapse load of the column was relatively small. 


3. For all aluminum alloy columns, the actual collapse load was found to 
be greater than the collapse load as predicted by the theory using the segment 
of a circle approximation but less than that using the cosine curve approxi- 
mation. Therefore, the circular segment is recommended for use in the de- 
sign of columns made from this and similar materials. 


4. For the six rectangular SAE 1020 steel columns subjected to sustained 
loads, the theoretical collapse loads were computed using the yield point 
stress as obtained from conventional compression and tension tests. In all 
cases the experimental collapse loads fell below the theoretical collapse 
load based on the segment of a circle approximation. Since this approxima- 
tion gives a conservative estimate of the collapse load, it can be concluded 
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that the yield point stress in members subjected to sustained loads is re- 
duced from 5 to 10 percent below the value obtained by standard testing 
procedures. 


5. Three rectangular steel columns were loaded through the inelastic 
range in time intervals totaling 6, 20, and 82 sec. For the latter two times, 
the rate of straining was comparable to that used in conventional tension and 
compression tests and the collapse loads were found to lie between the values 
computed using the two theoretical approximations. For the test completed in 
6 sec., a higher collapse load was obtained and it can be concluded that the 
yield point stress was raised from 5 to 10 percent above the value obtained 
from standard testing procedures. 
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Vill. APPENDIX A 


Proof that the Deflection of an Eccentrically-Loaded Column Varies 
Linearly Along the Straight-Line Portion of a Constant Depth of Yielding In- 
teraction Curve—In Fig. 16a is shown a member subjected to a compression 
load P and a bending moment M. Yielding has penetrated to a depth equal to 
ay on the compression side of the member while the strain on the tension 
side of the member is equal to the yield strain €,9. Thus the depth a of the 
cross-section is elastic. This condition of deformation corresponds to point 
A on the interaction curve shown in Fig. 16b. 

The stress distribution resulting from the load P and moment M is 
S-D-A-F-G-H. Instead of working with this stress distribution it is more 
convenient to represent the stress distribution by the rectangular distribution 
D-C-J-H plus the triangular distribution F-G-J minus the triangular distribu- 
tion F-C-A. The forces represented by these three distributions are 


1 
Fig = % 4K + 
1 


where Aj and Ag are the portions of the area that are inelastic and elastic, 
respectively, and K is the ratio of aj toa. The bending moment Mag is equal 
to the moment of these three forces about the centroidal axis. 


or 
M, = CE + 


where C is a constant for a given depth of yielding. The curvature of the 
member corresponding to the stress distributions S-D-A-F-G-H is given by 
Eq. 10 as 


1__ €e1* €e2 €e1 * €e2 
R kh a 


This is the curvature for the member corresponding to point A at the end of 
the straight line portion of the interaction curve (see Fig. 16b). For any 
point B on the straight line the stress on the tension side of the member is 
given by D-B and is less than 0e9. The resulting stress distribution for the 
same depth of yielding is U-D-B-F-T-H. This distribution is represented by 
the same rectangular distribution and two triangular distributions F-T-J and 
F-B-C. The moment arms m and n remain constant and, hence, the following 
expressions can be written: 
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The curvature corresponding to point B is given by the geometry of deforma- 
tion as (see Eq. 10) 


Comparison of the radii of curvature corresponding to points A and B shows 
that 


But it has been shown that 


and 


Therefore, 


Since the radius of curvature and the deflection are proportional, the de- 
flection corresponding to any point B on the straight line portion of the inter- 
action curve is given by the following relation: 
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THE PRESSURE LINE AND THE INELASTIC BUCKLING OF COLUMNS* 


Frank Baron* M. ASCE and Harold S. Davis** A.M. ASCE 
(Proc. Paper 1424) 


SYNOPSIS 


The pressure line concept developed for elastic arches and rigid 
frames,(1,2) and for the inelastic deformations of structural elements ds in 
length,(3) is applicable to studies of the inelastic behavior of columns and of 
rigid frames. The concept is extended herein for the inelastic behavior of 

such structures and for eccentrically and laterally loaded columns. In the 

same way as for structural elements, ds in length, the inelastic behavior of 
such columns can be stated in terms of the elastic behavior modified to take 
into account the effects of plasticity. 


INTRODUCTION 


The concept of a pressure line for a group of loads on an elastic structure 
is not new. For readers not familiar with the concept, a pressure line is for 
the present defined to be a string polygon that shows the successive resultants 
of the reactions and loads acting on a structure. The concept apparently has 
its origin in studies of masonry construction(4) and in the development of 
graphical methods of analysis. It has been used in studies of walls, 
,buttresses, and chimneys. It has also been used in studies of structures that 
lie in a plane and consist of single closed circuits. Examples of such struc- 
tures are arches and rigid frames. For these structures, the concept has 
ordinarily been used for elastic behavior and for loads that lie in the planes 
of the structures. 

The concept was recently extended(5) to structures that lie in a plane, are 
loaded normal to their plane, and are continuous between two supports. It 
was also extended(6) to structures that are curved in space, are loaded in any 


Note: Discussion open until March 1, 1958. Paper 1424 is part of the copyrighted 
Journal of the Engineering Mechanics Division of the American Society of Civil 
Engineers, Vol. 83, No. EM 4, October, 1957. 
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direction, and are continuous between two supports. In each case, however, 
elastic behavior only was considered. 

More recently, a procedure based on the pressure line concept was 
developed(7) for quickly estimating the effects of inelasticity on the behavior 
of structural elements, ds in length, subjected to axial and flexural loads. 
The procedure consists of obtaining an initial estimate by means of an ele- 
mentary theory of mechanics and adjusting the estimate to fit the conditions 
of a theory of plasticity. In this way, the moment-angle change relationship 
can readily be obtained for any value of axial thrust and moment. The proce- 
dure was demonstrated pictorially. It can also be used in conducting alge- 
braic studies. The latter is demonstrated herein for several kinds of stress- 
strain curves. ; 

The concept of the pressure line can also be used in studies of the inelastic 
behavior of arches and rigid frames and the inelastic buckling of eccentrically 
and laterally loaded columns. This is demonstrated herein in two ways; 
namely, pictorially for arches and frames and algebraically for columns. 


The Pressure Line for Arches and Frames 


This discussion is restricted to arches and frames that lie in a plane and 
are loaded in their plane. It is assumed that the loads are increased from 
zero to their final values. In addition, it is assumed that the effects of shear 
and residual strains can be neglected and that the dimensions of the unde- 
flected structure can be considered in determining the pressure line. It is 
further assumed that localized buckling, lateral-torsional buckling, or any 
other kind of buckling does not occur. Failure, if developed, is assumed to 
be caused by a sufficient number of hinges being formed to produce a mecha- 
nism in the plane of the arch. 

In 1879, Winkler of Berlin stated(8) in the “Zeitschrift des Architekten und 
Ingerieur Vereins zu Hannover,” page 199, a theorem concerning the position 
of the pressure line, or line of resistance, for single-span elastic arches. He 
stated, “For an arch ring of constant cross section that line of resistance is 
approximately the true one which lies nearest to the axis of the arch ring, as 
determined by the method of least squares.” This theorem was drawn from 
conclusions based on studies of the voussoir masonry arch and is only ap- 
proximately correct for the elastic theory of solid arches. 

The above theorem was extended by Hardy Cross, Hon. Mem. ASCE, in 
lectures at the University of Illinois and at Yale University for loads lying 
in the plane of single-span elastic arches and closed rings.(9) He restated 
the theorem in three forms as follows: 


1. “The pressure line for a group of loads lying in the plane of an arch is 
that string polygon for the loads which most nearly fits or hugs the axis of 
the arch.” (“Fit” in this statement is to be judged by eye.) 


2. “The pressure line for a group of loads is that string polygon which 
makes the moment areas balance or vanish.” (For vertical loads, the mo- 
ments in the arch are proportional to the vertical intercepts between the 
string polygon and the axis of the arch. For horizontal loads, the moments 
are proportional to horizontal intercepts. In generai the moment at any point 
in a loaded arch is equal to the force represented by the corresponding string 
of the pressure line times the perpendicular distance from the point to the 
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string. “Balance” in this statement means that the total positive area of the 
moment diagram equals the total negative area and that the centroid of the 
total positive area coincides with the centroid of the total negative area. This 
statement is for fixed-end arches of constant cross-section and can easily be 
modified for arches with hinges.) 


3. “The pressure line for a group of loads is that string polygon which 
makes the M/EI diagram balance or vanish.” (This statement is the same 
as the preceding statement except that provisions now are made for varia- 
tions in I along the axis of the arch.) 

For obtaining quick estimates of the effects of loads on an arch, Statement 
1 is the most useful of the above statements. Statements 2 and 3 successively 
modify the definition of “fit” and lead to a greater precision in computed ef- 
fects. Statement 3 fulfills the same requirements of statics and of geometry 
as are observed in the column analogy, the hydrostatic analogy, the neutral 
point method, and other formalized statements of the theory of the elastic 
arch. Statics is satisfied by specifying that the pressure line for a group of 
loads is a string polygon for the loads. All statical relationships in the arch - 
such as shears, thrusts, and moments - can be read directly from the string 
polygon. A force polygon for the loads is not needed in constructing the string 
polygon nor in obtaining the scale of moments in the arch. For example, the 
diagram of moments in an arch for a group of vertical loads on the arch is 
drawn to the same scale as that of the moments in a cantilevered beam or of 
the moments in a simple beam supporting the same loads and having the same 
span length as the arch. Also the intercept (parallel to a load and at a unit 
perpendicular distance away from the load) measured between the two rays of 
a string polygon intersecting at the load is proportional to the load and to the 
moment of the load about an axis at the intercept. The latter statement is 
particularly useful in obtaining all shears and thrusts in a loaded arch. 

The requirements of geometry for fixed-end arches loaded in the elastic 
range are fulfilled in Statement 3 by specifying that the M/EI diagram must 
balance. The requirements, stated in a more formal way, are as follows: 


(1a) 


(1b) 


(1c) 


For the inelastic behavior of arches and rigid frames, Statement 3 needs 
to be slightly modified to take into account the effects of non-linear stress- 
strain characteristics of materials. The modification is as follows: 


4. The pressure line for a group of loads on an arch is that string polygon 
for the loads which makes the M/E I diagram balance or vanish. In this 


statement, Ep is a modified modulus of elasticity and for an element one unit 
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in length is a measure of stiffness to relative rotation of the adjacent cross- 
sections. The value of E}, at a section depends on the shape of the cross- 
section, the stress-strain relationship, the history of loading, and on the 
magnitudes of the thrust and moment. For ultimate loads provisions can be 
made for the development of “plastic hinges.” 

For a specific case, the values or E}, can be obtained by means of the 
pressure line concept developed for the inelastic bending of structural 
elements.(10) This is illustrated in Fig. 1 for a rectangular cross-section 
subjected to given values of P and M. As in the usual theory of plasticity,(11) 
the distribution of strains across the section is assumed to be planar. In 
this paper, a history of loading is assumed that is statical. The strain in 
each fiber is assumed to be increased from zero to its final tensile or com- 
pressive value. Other histories of deformations are possible and are dis- 
cussed elsewhere.(12) The pressure line concept for other histories of defor- 
mations is illustrated in Ref. 3. In this paper, the stress-strain relationships 
for the fibers of a member subjected to a simultaneous axial and flexural load 
are assumed to be the same as those obtained from simple tension and com- 
pression tests. In Fig. 1, the stress-strain relationship of the material is 
represented by curve AOB and is assumed to be the same for each longitudinal 
fiber. As previously stated, the effects of a shear, Poisson’s ratio, residual 
strains, and local instability of projecting elements of the cross-section are 
herein neglected. In addition, the effects of lateral-torsional instability are 
not considered. 

In Fig. 1 the stress as given by the usual theory of plasticity is interpreted 
to be equal to the stress computed by the elementary theory of mechanics 
(So = Sq + Sp; where Sg = P/A and Sp = My-/I) plus a correction stress to ac- 
count for the non-linear stress-strain characteristics of the material. That 
is, 

SotSc (2) 


The ordinates bounded by the stress-strain curve AOB are in accordance with 
the theory of plasticity whereas those bounded by the straight line CD are in 
accordance with the elementary theory. The shaded diagram of Fig. 1 repre- 
sents the correction stresses as defined in Eq. 2. The correction stresses 
indicated in the figure must balance or vanish. “Balance” in this case means 
(1) the sum of the positive correction forces must equal the sum of the nega- 
tive correction forces, and (2) the centroid of the positive correction forces 
must coincide with the centroid of the negative correction forces. A correc- 
tion force is equal to a correction stress times its corresponding differential 
area. A pictorial procedure for effecting such balance has been presented 
elsewhere(13) and will not be described herein. The procedure is similar to 
that used in sketching the pressure line for fixed-end beams. 


Two modified moduli of elasticity are represented in Fig. 1. These are as 
follows: 


Eq=Sa/ég (3a) 
Eb =Asp/Ae (3b) 


where Sq = P/A, ASb = Mh/I, €q is the strain of the fibers along the cen- 
troidal axis, and Aé is the difference in the strains of the extreme fibers. 


4 
| 
: 
| 
‘| \ 
| 
4 
4 


ASCE 


BARON - DAVIS 1424-5 


The value of Ej, in this case is the slope of the line CD shown in the figure. 

The revised statement of the pressure line (Statement 4) used jointly with 
the above procedure for obtaining values of Ep is useful for quickly estimat- 
ing the inelastic behavior of structures such as arches and rigid frames. The 


suggested procedure for making these estimates consists of the following 
steps: 


1) Draw a string polygon for a given set of loads which you think is the 
best “fit” to the axis of the given arch. Obtain scale on the magnitudes of the 
moments as interpreted through the string polygon. 


2) Check “fit” by inspecting whether the moment diagram balances. Re- 
vise, if necessary, the string polygon to obtain a more satisfactory balance. 


3) Draw the M/EI diagram and check the balance of this diagram. Revise 
the string polygon until a satisfactory balance of the M/EI diagram is obtained. 
(The revision is aided by noting that a section with a small value of I tends 
to draw the pressure line toward itself whereas a section with a large value of 
I tends to push the pressure line away. In other words, a section with a small 
value of I tends to attract less moment whereas a section with a large value 
of I tends to attract more moment). Now by means of the string polygon, ob- 
tain all shears, thrusts, and moments in the arch. (Note that the results at 
this stage are in accordance with the elastic theory based on the original 
dimensions of the structure.) 


4) For the distribution of moments and corresponding thrusts obtained 
above, estimate the values of Ep: (If convenient, draw diagrams relating 
moments with angle-changes for different values of thrust. See Fig. 5a) 


5) Draw the M/EpI diagram and check the possible unbalance. If neces- 
sary, revise the string polygon and the values of Ep until the M/EpI diagram 
is considered sufficiently balanced. (The revision here is aided by noting that 
a section with a reduced modulus of elasticity tends to draw the pressure line 
toward itself.) Now interpret all shears, thrusts, and moments by means of 
the string polygon. 

Step 3 is not necessary in studies restricted to inelastic behavior. How- 
ever, it is useful in developing judgment as to the final position of the pres- 
sure line. It also is useful in comparing the results of the elastic theory with 
those of the inelastic theory. This in itself suggests that interaction dia- 
grams, as for columns with axial and flexural loads, can be developed for the 
elastic and inelastic behavior of arches and rigid frames. 

In the above procedure, the distributions of angle-changes per unit of 
length are assumed to be known for both theories. Consequently, estimates 
can be made of the corresponding deflections. If desired, the above procedure 
can be extended to the theory of the flexible arch. The pressure line in that 
case must fit the axis of the deflected arch. For these studies, larger scales 
usually are required for the sketches of the arch and more accuracy in the 
computations of fit. The latter computations however become sensitive to 
changes in the modified moduli of elasticity resulting from changes in mo- 
ments and deflections of the arch. This is particularly true when total yield- 
ing is approached. For these cases, the sketching procedure becomes inade- 
quate and a numerical procedure of successive approximations, similar to 
that of Vianello,(14,15) is usually required. 
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A numerical example of a rigid frame with fixed ends and a concentrated 
load, W, at mid-span is shown in Fig. 2. The dimensions and properties of 
the frame are given in the Fig. A pressure line also is drawn in the Fig. for 
each of three ranges of loading. In Range 1, the load is 105 kips and produces 
a maximum stress of 33 ksi in an extreme fiber of the section at mid-span. 
The distributions of moments and angle-changes are shown and the require- 
ments of geometry are checked. In Range 2, the load is 171 kips and produces 
a plastic “hinge” at the mid-span of the frame. The maximum moment, My, 
that the given cross section can resist is 594 K ft. A plastic zone extends 
over a portion of the horizontal member as shown in the Fig. The distribution 
of angle-changes along this portion was obtained from a diagram relating mo- 
ments and angle-changes as shown in Fig. 5a. For the rise to span ratio con- 
sidered, the influence of the thrust on the moment-angle change relationship 
is negligible and consequently is neglected. The requirements of geometry 
for Range 2 are checked as for Range 1 except that provisions now are made 
for the effects of plasticity. Compare the pressure line for Range 2 with that 
for Range 1 and observe that the plastic zone attracts the pressure line to- 
ward itself whereas the sections at the abutments and the knees push the 
pressure line away. 

In Range 3, the load is 216 kips and produces yielding at a sufficient num- 
ber of sections to cause collapse. This load was obtained by means of the 
pressure line as well as by means of the “mechanism” criterion discussed 
elsewhere.(16) For the pressure line drawn in the Fig., observe the extent 
to which it hugs the axis of the frame. As for Range 2, the influence of the 
thrust on the maximum moment that a section can take is negligible and 
consequently was discounted. 

For small rise to span ratios, the influence of the thrust on the moment- 
angle change relationship can be appreciable and merits consideration. In 
addition, the influence of deflections on moments is usually important for 
such ratios and needs to be taken into account. 

In the above example, no consideration was given to lateral-torsional 
instability. It is assumed that the frame is sufficiently braced for such in- 
stability not to occur. Obviously, if such instability can occur, the limiting 
load should be less than that obtained as the moment capacities of the “plastic 
hinges” could not be developed. 


Inelastic Buckling of Columns 


The inelastic behavior of a column can be studied in the same way as above 
except that the pressure line for a group of loads must now be in agreement 
with the final deflections of the column. In addition, the curvatures of the 
column, proportional to Ag = M/EpI, must be in agreement with the values 
of Ep and with the deflections of the column. As indicated in Fig. 1, the 
value of Ep at a section depends on the shape of the cross-section, the 
stress-strain relationship, and on the magnitudes of the thrust and moment. 
The value of E}, also depends on the history of loading. The history of load- 
ing and the mode of failure are of such importance in column theory that cer- 
tain assumptions made in the preceding discussion are herein repeated. The 
strain in each fiber is assumed to gradually increase from zero to its final 
tensile or compressive value. Other histories of deformations are possible, 
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including residual strains, and are discussed elsewhere.(17) Local instability 
and lateral-torsional instability are not considered. The latter type of 
instability can be important and is frequently obtained in laboratory tests.(18) 
In such tests, no resistance is frequently offered to a twist at the ends or to a 
yaw and a twist at the mid-portion of a column. In actual practice, such 
resistance is frequently offered by end connections and by walls or bracing 
along the length of the column. The succeeding discussion is based on the 
assumption that the failure of a column is due to excessive bending about the 
axes of moments. The discussion illustrates the use of the pressure line con- 
cept for such columns. The development is analytical rather than pictorial 
or numerical. 

Fig. 3 shows some typical pressure lines (curves of moments) for various 
conditions of loading on columns. It is assumed for the present that the loads 
are less than the critical loads at which instability occurs. In each diagram, 
curve e represents the deflected position of the column as if the behavior of 
the column was completely elastic. A generalized statement concerning the 
behavior of such columns was given by the late H. M. Westergaard.(19) The 
statement was given in the form of a Fourier series and considerable liberty 
is taken by the writers in restating it as follows: The elastic behavior (later- 
al deflections and moments) of a structural member subjected to axial and 
lateral loads is equal to its behavior as a laterally loaded member multiplied 
by a magnification factor. The magnification factor depends on the ratio of 
the axial load to the Euler load for an idealized straight column. 

In Fig. 3, a comparison also is shown between the elastic and the inelastic 
behavior of each column. Curve i in each diagram represents the deflections 
of the column in accordance with a theory of inelasticity. The inelastic be- 
havior of the column can be stated in the same way as the elastic behavior 
modified to take into account the effects of plasticity.(20) This is illustrated 
for the case of a column as shown in Fig. 3c. The column is assumed to have 
a symmetrical distribution of lateral loads along its length and to have a load 
P of equal eccentricity at each end. The initial eccentricity is designated by 
Yo: It is assumed that the column has a symmetrical curve of deflections, 
that bending occurs about axes normal to the plane of the bent column, and 
that the cross-section of the column has an axis of symmetry lying in this 
plane. The final deflections at the center of the column for conditions of 
elastic and inelastic behavior respectively are designated y, and yj. The 
effects of residual strains, lateral and torsional instability, and local buckling 
are not considered. In addition, uncertainties in the properties of materials 
are not a part of this discussion. 

In the succeeding discussion, the elastic and inelastic behavior of the 
column will be paralleled and will respectively be associated with the sub- 
scripts e and i. The final deflections at the center of the column are defined 
by 


A. Elastic B. Inelastic 


ye = Yo tYae Yi = Yot yai (4) 


where y,, and y,; are the additional deflections caused by the loads. The 


pressure line for the group of loads must coincide with the resultant of the 
stress distribution across the section. By definition, the moments at the 
center of the column are given by 
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A. Elastic B. Inelastic 


Me=M,+P(yo+yae) Mi=Mi+P(yotygi) (5) 


where My is the moment due to the lateral loads. Eqs. 5 can be rewritten as 


A. Elastic B. Inelastic 


Me=Momt+Mge Mj= Mom+Mai (6a) 


Mae= Pyae Mgi = Pyai (6b) 


Mom =[!+M, /Me]Mo (6c) 


and where Mb) = Pyo for elastic and inelastic behavior. It is observed that 
Mom is the moment at the center of the column with respect to the original 
dimensions of the column and is a constant for a given case of loading. 

The additional deflections caused by the loads can be written as follows: 


A. Elastic 


Mom=[!+M_/Mo|Mo 


B. Inelastic 
oe 2 


where ky and kj are coefficients depending on the distributions of the angle- 
changes (or curvatures) along the length of the column and A, and Ag; are 


the angle changes per unit of length at the center of the column. The angle- 
changes per unit of length are related to the moments as defined by 


A. Elastic B. Inelastic 


Mi (8) 
Abe= ET 


where E is the tangent modulus at the origin of the stress-strain curve and 
E), is a reduced modulus of elasticity. 

By means of the above relationships and definition of the Euler load, 
PR = n2EI/L2, the deflections at the center of the column are as follows: 


A. Elastic 


B. Inelastic 
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| 
ye* Yor! Yom yi= Yor Pe |Yom (9b) 
ep ci 
where 
Ke E 
Ce“ TRE cj= (90) 


and where Yom = (1 + M,/M,) Yo for each behavior. Consequently, M,., 
= Pyom- In the same way as for deflections, the moments at the center of 
the column are given by 


A. Elastic B. Inelastic a 
| 10 
Mae= Pe. \Mom q 
Ce Ci 
M Ce “5 M Cj Pe (10b) 
P om Mj = wi Mom 


In Eqs. 9a to 10b, the terms in the brackets can respectively be interpreted 
as modification factors of yoy, & Mom. A comparison of the moments de- 
fined by Eqs. 10a and 10b yields 


A. Elastic B. Inelastic X 

Me Pe Mj (11) 

Mage e Magi 4 


It is of interest to note that the form of Eqs. 4 to 11 are the same for elastic 
and inelastic behavior except for the values of cg and cj. A further compari- 


son of the moments defined by Eqs. 10a and 10b yields the following relation- 


ship: 
Mi /Me _| ci ki Ep 
Mai/ Mae | Ke E 


Thus, in a general way the inelastic behavior is the elastic behavior modified 
by the terms in the brackets of Eq. 12. 
Values of ce for elastic columns have been reported by H. M. 


Westegaard(21) for various conditions of axial and lateral loads. Fora 2 a 
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symmetrical condition of load, the value of c, is not sensitive to the type of 
loading and is approximately equal to 1.0; that is, k, = m2, The values of Cj 
in the above Eqs. depend on the deflected shape of the column and on the 
P-M- A¢ relationships for a given cross-section and stress-strain charac- 
teristics of the column. 

The values of Ep/E in the expression for cj can be obtained by means of the 
pressure line concept developed for the inelastic deformations of structural 
elements. For example, consider an element with a rectangular cross- 
section and an idealized stress-strain curve as shown in Fig. 4. The stress- 
strain curve consists of three straight lines defining material properties 
which are the same for tension as for compression. The slope of the line 
through the origin is designated by E and that of the other lines by Et, 
(where n = 1 - E;}/E). For reference, the coordinates of the intersections of 
these lines are defined by Sy, fy and “By, - €y- Other shapes of cross- 
sections and stress-strain relationships can be considered. These are dealt 
with elsewhere.(22) 

In Fig. 4, several ranges of loading (1 to 4) are defined. In Range 1, the 
behavior of a column is elastic, whereas in the other Ranges the behavior is 
inelastic. In Range 2, only the compressive stresses on the concave side of a 
column exceed Sy- In Range 3, the compressive stresses on the concave 


side exceed Sy and the tensile stresses on the convex side exceed -Sy. 


In Range 4, all stresses are compressive and exceed Sy- The shaded 
diagrams in the Figure represent the correction stresses to the elemen- 


tary theory of mechanics as defined by Eq. 3. The slope of line CD for each 
range of loading represents the corresponding value of the modified modulus 
of elasticity, E,. The expressions for Ep can be obtained by observing that 
the shaded diagrams must balance. (23) They can also be obtained by ob- 
serving that the correction stresses at the extreme fibers of the cross- 
section are numerically the same as the fixed-end moments of a fictitious 
beam having a span-length of t and having concentrated loads equal to E-E; 
at the corresponding positions of the abrupt changes in the stress-strain curve. 
This is illustrated in Fig. 4. 

In the succeeding discussion, particular consideration is given to a 
material having a well defined yield point; that is, n = 1 -E;/E = 1. Signifi- 
cant equations associated with the generalized stress-strain curve are sum- 


marized in Fig. 4. When n = 1, the equations for Ep/E and Ad/A by are as 
follows: 


Range 2: (n = 1) Range 3: (n = 1) 


I-PP I-P/Py) E 
Eb_M [5(I- Eb M 
E My 4(I-P/A)3 y 

—4(I-P/PyP Ad ese 
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where the subscripts i previously associated with Ad; and M; have been 
dropped for convenience in writing. Families of curves,(24) based upon the 
relationships defined by these Eqs. are shown in Fig. 5. 

In Fig. 5, curves also are given which relate the maximum values of M/ 
for WF sections when the values of P/P, are kept constant. More complete 


curves and corresponding expressions for P - M - Ad relationships of WF 
and other rolled sections are given in Reference 20 and elsewhere.(25,26) 

In Figure 6, curves are given(27) which relate P/Py, M/My, and A¢/A dy for 
rectangular cross-sections when E,E = 0.1. The influence of the stress- 


strain characteristics of a material on the resulting P - M - A¢@ relationships 
can be seen by comparing Fig. 5a with Fig. 6b. 


Critical Loads 


Ros(28) defined a critical load to be obtained when the external moment at 
any section of a column increases at a rate equal to the rate of increase of the 
internal resisting moment offered at the same section. This criterion stated 
in the form of an equation is as follows: 


d(M)| |d(m) 
ext Had) int. 


It is used herein(29,30) to develop significant expressions associated with 
critical loads. Since the external moment (See Fig. 3) at a section can be 
written as 


Mext, = M_+P(yo+ Yai) (15) 


its rate of increase becomes 


; y a This assumes that kj is a constant during the increase. The critical load Z 
: then becomes 


This permits a more direct means for determining critical values of slope 
along P - M -A®@ curves than semi-graphical methods used in the past. 
Proceeding as previously, we obtain by means of Eqs. 13, the following 


Eqs. for [aca 4 ¢) | int? D2mely; 


eee 
ag 
d(M)| _ d(Pyqi) __d(Pyai)_, 2 (16) 
d(a 6) d(ad) d K Kj 
4 
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d(M) 
Substituting these Eqs. into Eq. 17 and noting that a 


Agy Ep 


the following relationships associated with critical loads are obtained: 


Range 2: (n = 1) Range 3: (n = 1) 


(19a) 


. (19b) 


(19c) 


Eb 3 
E cr Pe kj 
The relationships defined by these Eqs. are represented graphically in Fig. 7. 


The relationships shown in Fig. 7a are of particular interest as they consist 
of straight lines for P..../Pg equal to a constant. 


Interaction diagrams, as suggested by Shanley,(31) can now be determined 
by means of the preceding equations and 


Substitution into this Eq. of the relationships defined by Eqs. 7, 19a, and 19c 
yields the following results: 


Range 2: (n = 1) Range 3: (n = 1) 


These Eqs. can be rewritten in terms of L/r, since 


*See Fig. 4 for the corresponding equations when n is not equal to 1.0. 
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for each range of loading. The resulting expressions define interaction dia- 
grams as drawn in Fig. 8. Values of k, = 72 and Sy/E = 33/30,000 were as- 
sumed in plotting this Figure. 

The eccentricity ratio, Yom/h2r2, associated with the value of a critical 


axial load on a laterally loaded column can readily be determined in Fig. 8 by 
observing that 


Yom h_ M P 


Thus, the inverse value of the slope of a line drawn from the origin of the 
Fig. is equal to the eccentricity ratio associated with the corresponding 
values of Por/Py and Mom/My. The intersection of such a line in Fig. 8 
with a curve for L/r equal to a constant defines the corresponding values of 
Por/Py and Mom/My of the loaded column. This same concept has been used 
in constructing interaction diagrams for the flexural stresses and deflections 
of flexible arch ribs(32) and of the stiffening girders of suspensions 
bridges.(33) Similar constructions can be made in several of the other Figs., 
as for example in Fig. 5a where 


Ex _M (24) 


The values of kj in the preceding equations are not sensitive to reasonable 
assumptions concerning the distributions of angle changes along a column.(34) 
This is particularly true for straight columns loaded eccentrically at their 
ends and for columns resisting axial and lateral loads. The values of kj for 
triangular, parabolic, sinesoidal, and rectangular distributions of angle- 
changes along a pin-ended column respectively are 12, 9.6, 9.87, and 8. For 
a straight column loaded eccentrically at its ends, the value of k; is usually 
between 8 and 7r2. However, for rectangular columns which are initially 
bent, studies reported elsewhere(35) show that the values of kj may vary 
from around 74 to 16. For these cases, values of k; equal to 7%, 12 and 16 
correspond to columns having slenderness ratios of about 100, 40 and 15, 
respectively. 

For k; equal to 72, Eqs. 21 yield the same results as those interpreted 
by means of the Eqs. developed by Jezek.(36) The Eqs. herein, however, are 
more general and permit allowances to be made for variations in k,. It also 
is observed that the methods and concepts presented herein are general and 
can be used to determine P - M - A¢ relationships and interaction diagrams 
for coluntns having other shapes of cross-sections and stress-strain curves 
than those presented in this paper. 

Again, no consideration has been given in the preceding discussion to the 
possibility of lateral-torsional instability. Consequently the expressions and 
diagrams that are given are only valid when failure occurs due to excessive 
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bending about the axes of moments. It is recognized that if lateral-torsional 
instability can occur the critical loads for eccentrically and laterally loaded 
columns can be less than those given herein. For example see Fig. 8. 


CONCLUSIONS 


The concept of the pressure line is useful for estimating the inelastic be- 
havior of structural frameworks. It is equally useful for making analytical 
studies of the inelastic behavior of columns. 

In addition, the inelastic behavior of columns for the cases studied can be 
stated in terms of the elastic behavior modified to take into account the ef- 
fects of plasticity. 
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(b)Eccentric Loads (c)Eccentric & Lateral 


(a) Initially Crooked 
Loads 


LOADS AND INITIAL DIMENSIONS 


INELASTIC 


Defiections Moments 


ELASTIC 


Deflections Moments 
Ye * Yo t+ Yoe Me = M.+Ply¥o+yoe) Yi *YotYoi Mi = +P(yot yoi) 
Me *Mom +Moe Mi * Mom 


Moe iMo=P¥o Yai * 


FINAL DEFLECTIONS AND MOMENTS 


FIG.3_ ECCENTRICALLY AND LATERALLY LOADED COLUMNS 
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FIG. 4_ RELATIONSHIPS FOR GENERALIZED STRESS-STRAIN CURVE 
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FIG.6_RELATIONSHIPS BETWEEN PM,AQ AND E, WHEN n=0.9 
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FIG.7- RELATIONSHIPS BETWEEN AND FR AT VALUES OF CRITICAL LOADS 
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FIG. 8-INTERACTION DIAGRAMS FOR ECCENTRICALLY AND LATERALLY 
LOADED COLUMNS 


Lo — 
P 
i 
os 
Por 3.0 
Py 08 04 0.2 
04 4.0 
0.2 
| 04 O6 1.0 1.2 .4 L6 
(a) 
| Yorn * 05 Lo 5 
L/r 
20 
Pp 
5-€ Curve al 
60 
Py 80 
0.4 40 
a2 
0 
4 
| 


i 
| 
i: 
| 
i 
| 
| 
1 
| 
| 
oy 


95 


83 


The) 
— 
a 
GS 
‘on 


4 


